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I.  ABSTRACT 
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group  \$ub-qroup  -^Liquefacti on ,  Porewater  Pressure,  Dynamic  Testing,  Blast 

_  f  Loading,  Soil  Mechanics,  Laboratory  Testing,  Field  Testing, 

1  vi  Material  Modeling,  Explosive  Loading,  Saturated  Sand.^Jg5i 

'  ~'~‘J0ur  field,  laboratory  and  theoreticalresearch^ndicate,  that  the 

destruction  potential  of  an  explosion  may  be  greatly  magnified  if 
detonated  in  water  saturated  granular  soils.  While  blast-induced 
liquefaction  may '.not  necessarily  damage  a  facility  structurally,  it  may 
render  it  unusable. \  Empirical  models  are  given  that  can  be  used  to 
estimate  liquefaction  potential  as  a  function  of  density,  effective  stress 
and  applied  compressive  strain.  One  of  the  models  uses  an  empirical 
scaling  law  for  explosive  loadings  to  predict  the  extent  of  porewater 
pressure  increases  in,  the  field  from  buried,  contained  charges  In 
saturated  soils.  A  numerical  analysis  that  considers  the  saturated  soil 
as  a  two-phase  medium  is  presented.  The  analysis  accounts  for  the 
nonlinear,  inelastic  behavior  of  the  soil  skeleton  and  has  shown  that 
liquefaction  is  dependent  upon  the  constrained  modulus  of  the  soil 
skeleton.  ^Results  agree  with  the  experimental  observations.  (cont  ) 
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{Abstract  Continued) 


, The  results  of  our  study  Indicate  the  following, 

1.  Liquefaction  can  be  induced  by  single  and  multiple  blasts. 

2.  Liquefaction  can  be  induced  at  distances  much  greater  than 
those  associated  with  structural  damage. 

3.  Long  term  increasas  in  residual  porewater  pressures  can  be 
induced  by  compressive  shock  wave  loadings  when  the  peak 
particle  velocity  exceeds  0.075  m/s. 

4.  Liquefaction  can  be  induced  in  loose  saturated  sand  by  a 
single  compressive  shock  wave  when  the  peak  particle  velocity 
exceeds  0.75  m/s. 

5.  Soils  at  higher  initial  effective  stress  and  higher  initial 
relative  density  require  more  energy  to  produce  liquefaction. 

6.  Destruction  potential  of  an  explosive  charge  may  be  greatly 
magnified  if  detonated  in  water-saturated  soils. 

7.  Liquefaction  occurs  because  of  compressive  strain  induced  by 
the  compression  stress  wave,  but  liquefaction  occurs  after  the 
stress  wave  passes. 

The  possibility  exists  that  an  explosive  detonated  in  a  soil 
having  a  high  liquefaction  potential  could  cause  damage  disproportionate 
to  the  energy  released.  Documented  occurrence  of  blast-induced 
liquefaction  is  available  in  the  open  literature.  Although  considerable 
work  remains  to  be  done  in  projecting  this  information  into  a 
comprehensive  method  of  predicting  liquefaction  for  actual  or  hypothetical 
blasts,  the  data  indicate  that  residual  porewater  pressure  Increases 
should  not  occur  in  soils  subject  to  strains  of  less  than  0.*005  percent. 
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I.  INTRODUCTION 


A.  Blast  Induced  Liquefaction  Threat 

Blast- induced  liquefaction  of  water  saturated  soils  represents  a 
threat  to  both  military  and  civilian  structures,  liquefaction  can  lead  to 
catastrophic  consequences  including  landslides,  foundation  failure, 
floatation  of  buried  buoyant  structures,  ground  subsidence,  and  failure  of 
earthfill  dams.  Liquefaction  may  be  the  cause  of  the  unusually  broad, 
flat  crater  shapes  and  late-time  low  frequency  ground  motion  observed  in 
nuclear  (N£)  and  high  explosive  (HE)  tests  in  the  Pacific  and  Canada, 
respectively. 

At  the  International  Workshop  on  Blast-Induced  Liquefaction, 
McCracken  (1978)  stated  that,  ,lWe  within  the  United  States  Air  Force  now 
believe  that  blast-induced  soil  liquefaction  could  be  a  far  more  serious 
threat  to  both  civilian  and  military  targets  than  we  have  given  it  credit 
for  in  the  past."  The  Dutch  also  consider  liquefaction  a  serious  threat 
as  is  evident  by  Kok's  (1978)  statement  during  a  discussion  on  possible 
effects  from  the  detonation  of  a  500,000  kilogram  (500  kiloton)  nuclear 
weapon.  Kok  stated  that,  "everything  will  be  liquefied  ...  all  the 
structures  will  fall  to  pieces."  Rischbieter  (1977)  noted  that  the 
Netherlands  are  concerned  because  of  their  coastal  plains  lying  below  sea 
level  and  Switzerland  is  concerned  because  of  their  abundance  of  lake 
shore  deposits  consisting  of  post-glacial  alluvium.  In  Norway,  Kummeneje 
and  Eide  (1961)  showed  that  excess  porewater  pressure  and  liquefaction 
could  be  induced  by  blasting.  The  Russians  have  conducted  an  extensive 
series  of  field  explosive  tests  and  have  developed  empirical  methods  to 
predict  the  extent  of  blast-induced  liquefaction  (Lyakhov,  1961;  Ivanov, 
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1967;  and  Florin  and  Ivanov,  1961).  One  of  the  earlier  references  to 
liquefaction  was  in  a  paper  by  Terzaghi  (1956)  in  which  he  stated  that 
nearby  blasting  operations  caused  the  1935  failure  of  the  SWIR  II  dam  in 
Russia. 

B.  Definition  of  liquefaction 

In  1978,  The  American  Society  of  Civil  Engineers,  Committee  on 
Soil  Dynamics,  Geotechnical  Engineering  Division,  defined  liquefaction  as. 

The  act  or  process  of  transforming  any  substance  into  a  liquid. 
In  cohesionless  soils,  the  transformation  is  from  a  solid  state  to 
a  liquified  state  as  a  consequence  of  increased  pore  pressure  and 
reduced  effective  stress. 

Several  other  definitions  of  liquefaction  exist  in  the  literature. 
Basically,  they  all  state  that  liquefaction  is  a  condition  caused  by  an 
increase  in  porewater  pressure,  thus,  a  decrease  in  effective  stress  and  a 
loss  of  shear  strength  of  a  soil  mass. 

Liquefaction  can  be  caused  by  several  mechanisms  which  include  any 
activity  in  saturated,  cohesionless  soil  that  cause  the  soil  grains  to 
form  a  more  compact  structure.  This  compaction  results  in  an  increase  in 
porewater  pressure  which,  in  turn,  leads  to  a  decrease  in  effective  stress 
and  thus,  liquefaction.  Liquefaction  can  last  for  seconds,  minutes,  or 
hours  and  longer.  Excess  porewater  pressure  has  been  reported  to  last  up 
to  several  days.  Given  sufficient  time,  gravity  induced  failures  can 
occur. 

Liquefaction  effects  take  many  forms.  These  include  flow  failures 
of  slopes  or  earth  dams,  settlement  or  tipping  of  buildings  and  piers, 
collapse  of  retaining  walls,  lateral  spreading  of  inclined  ground,  and 
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defomation  of  the  ground  surface.  Water  spouts  and  sand  boils  typically 
accompany  liquefaction. 

C.  Objectives  of  the  Research 

The  primary  objective  of  our  investigation  was  to  systematically 
evaluate  the  behavior  of  saturated  granular  soils  subjected  to  shock  and 
explosive  loadings.  We  conducted  laboratory  and  field  experiments  to 
simulate  the  field  loading  of  a  soil  element  located  near  the  detonation 
point  of  an  explosive.  Intense  compressive  loadings  having  millisecond 
rise  times  to  peak  stress  occur  in  this  region.  The  soil's  porewater 
pressure  response  both  during  and  after  the  passage  of  the  stress  wave  was 
used  to  evaluate  the  liquefaction  potential  of  the  soil.  Experimental 
measurements  included  the  applied  loading  stress,  peak  particle  velocity 
and  the  porewater  pressure  response. 

A  large  number  of  parameters  have  been  observed  to  affect  the 
onset  of  liquefaction.  Some  are  associated  with  the  soil  while  others  are 
related  to  the  explosive  itself.  We  investigated  the  effect  of  variations 
in  the  initial  relative  density,  the  initial  effective  stress,  and  soil 
type  along  with  the  intensity  geometry  and  number  of  applied  loadings. 

Laboratory  shock  tests  were  conducted  on  two  quartz  beach  sands 
(coarse  and  fine),  a  granitic  river  sand,  a  granitic  river  gravel,  a  coral 
beach  sand,  and  a  clayey  silt.  Field  explosive  tests  were  conducted  on  a 
granitic  river  sand,  and  an  in-situ  granitic  river  sand-gravel  deposit. 
Our  analysis  of  the  data  include  an  evaluation  of  the  influence  of  several 
parameters  on  the  peak  and  long-term  porewater  pressure  response  in  the 
soil,  the  stress  wave  velocity,  and  the  peak  particle  velocity.  The 
results  of  these  analyses  are  used  to  define  liquefaction  threshold  limits 


and  develop  empirical  relationships  for  predicting  porewater  pressure 
increases  in  saturated,  cohesionless  soils.  Our  study  has  derived  and 
documented  several  important  relationships  between  soil  properties  and 
compressional  stress  wave  loading.  In  the  course  of  this  work,  we 
evaluated  existing  empirical  predictive  techniques  including  the 
relationships  between  scaled  charge  distance,  peak  particle  velocity,  and 
peak  porewater  pressure,  versus  porewater  pressure  response.  We  use 
theoretical  relationships  and  analytical  models  to  explain  the  observed 
behavior  of  our  field  and  laboratory  tests. 

0.  Organization  of  the  Report 

In  Chapter  II,  we  present  a  review  of  the  literature  on  blast- 
induced  liquefaction.  Chapters  III  and  IV  present  the  results  of 
laboratory  shock  tests  on  water  saturated  gravels,  sands  and  clayey  silts 
conducted  by  Veyera  (1985),  Hubert,  (1986),  Chouicha  (1987)  and  Bolton 
(1988).  Chapter  IV  also  presents  the  results  of  laboratory  high-pressure 
quasi -static  tests  on  sand,  clayey  silt  and  sand-clayey  silt  mixtures 
conducted  by  Bolton  (1988).  Chapters  V  and  VI  present  the  results  of 
field  explosive  tests  on  saturated  sands  conducted  by  Bretz  (1988),  Hassen 
(1988),  Jacobs  (1988),  Schure  (1988)  and  Allard  (1988).  Chapter  VII 
presents  empirical,  analytical  and  theoretical  analysis  of  shock  and 
explosive  induced  liquefaction  conducted  by  Veyera  (1985),  Awad  (1988), 
Hassen  (1988),  Bretz  (1988)  and  the  authors  of  this  report. 
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E.  Conclusions 

Our  research  indicates  that  the  destruction  potential  of  an 
explosion  may  be  greatly  magnified  if  detonated  in  water  saturated 
granular  soils.  While  blast-induced  liquefaction  may  not  necessarily 
damage  a  facility  structurally,  it  may  render  it  unusable.  Blast- induced 
liquefaction  can  cause  late  time  decreases  in  the  soil's  shear  strength 
that  produces  damage  disproportionate  to  the  amount  of  explosive  used  and 
ground  motions  inconsistent  with  previous  experience.  For  example,  recent 
re-examinations  of  the  events  at  the  Pacific  Proving  Grounds,  where 
nuclear  explosives  (ME)  were  detonated  in  the  1950's,  and  high  explosive 
(HE)  tests  conducted  at  Suffield,  Canada,  suggest  that  liquefaction  may  be 
the  primary  factor  causing  the  unusually  broad,  flat  crater  shapes. 

The  results  of  our  study  indicate  the  following. 

1.  Liquefaction  can  be  induced  by  single  and  multiple  explosive 
induced  compressive  wave  loadings. 

2.  Liquefaction  can  be  induced  at  distances  from  explosions  much 
greater  than  those  associated  with  structural  damage. 

3.  Fairly  long  term  increases  in  residual  porewater  pressures  can 
be  induced  by  compressive  shock  wave  loadings  when  the  peak 
particle  velocity  exceeds  0.075  m/s,  the  peak  porewater 
pressure  exceeds  250  kPa,  or  the  peak  strain  exceeds  0.005 
percent. 

4.  Liquefaction  can  be  induced  in  loose  saturated  sands  by  a 
single  compressive  shock  wave  when  the  peak  particle  velocity 
exceeds  0.75  m/s,  the  peak  porewater  pressure  exceeds  2,500 
kPa,  or  the  peak  strain  exceeds  0.05  percent. 
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5.  Soils  at  higher  initial  effective  stress  and  higher  initial 
relative  density  require  more  energy  to  produce  liquefaction. 

6.  Destruction  potential  of  an  explosive  charge  may  be  greatly 
magnified  if  detonated  in  water-saturated  soils. 

7.  Liquefaction  occurs  because  of  compressive  strain  induced  by 
the  compression  stress  wave,  but  liquefaction  occurs  after  the 
stress  wave  passes. 

8.  Liquefaction  occurs  because  loading-unloading  of  the  porewater 
is  elastic  but  the  soil  skeleton  is  not. 

An  explosive  detonated  in  a  soil  having  a  high  liquefaction 
potential  could  cause  damage  disproportionate  to  the  energy  released. 
Documented  occurrence  of  blast-induced  liquefaction  is  available  in  the 
open  literature.  Although  considerable  work  remains  to  be  done  in 
projecting  this  information  into  a  comprehensive  method  of  predicting 
liquefaction  for  actual  or  hypothetical  blasts,  the  data  indicate  that 
residual  porewater  pressure  increases  should  not  occur  in  soils  subject  to 
strains  of  less  than  0.005  percent.  Transient  and  quasi-static  tests 
indicate  that  residual  porewater  pressure  increases  and  liquefaction  are 
not  strain-rate  sensitive  for  sand  but  are  strain-rate  sensitive  for  silt. 
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II.  LITERATURE  REVIEW 

This  chapter  summarizes  the  current  state  of  knowledge  about  the 
nature  and  occurrence  of  blast-induced  liquefaction  from  field 
observations,  from  laboratory  tests,  and  from  small-scale  field  tests. 
Empirical,  analytical  and  theoretical  aspects  of  liquefaction  are 
reviewed. 


A.  Introduction  and  Definition  of  Liquefaction 

The  term  liquefaction  has  been  used  to  describe  the  state  of  a 
saturated,  cohesionless  soil.  The  manifestations  of  soil  in  a  liquefied 
state  are:  sand  boils,  flow  failures,  low  frequency  ground  oscillations, 
loss  of  bearing  capacity,  rise  of  buoyant  buried  structures  and  ground 
settlement.  Indirect  evidence  of  a  liquefied  soil  include  the  measurement 
of  excess  porewater  pressure  and  delayed  failures  (Committee  on  Earthquake 
Engineering,  1985). 

Several  definitions  have  b°en  proposed  for  liquefaction.  These 

include: 

Liquefaction-  “Denotes  a  condition  where  a  soil  will  undergo 
continued  deformation  at  a  low  residual  resistance  stress  or  with 
no  residual  resistance  due  to  the  build  up  of  high  pore-water 
pressures  which  the  effective  confining  pressure  to  a  very  low 
value;  pore-water  pressure  build  up  leading  to  true  liquefaction 
of  this  type  may  be  due  to  either  static  or  cyclic  stress 
applications"  (Seed,  1975). 

Initial  Liquefaction-  "Denotes  a  condition  where,  during  the 
course  of  cyclic  stress  application,  the  residual  pore-water 
pressure  on  completion  of  any  full  stress  cycle  becomes  equal  to 
the  applied  confining  pressure;  the  development  of  initial 
liquefaction  has  no  implication  concerning  the  magnitude  of  the 
deformations  which  the  soil  might  subsequently  undergo;  however  it 
defines  a  condition  which  is  a  useful  basis  for  assessing  various 
possible  forms  of  subsequent  soil  behavior"  (Seed,  1976). 

Liquefaction-  "The  act  or  process  of  transforming  any  substance 
into  a  liquid.  In  cohesionless  soil,  the  transformation  is  from  a 
solid  state  to  a  liquefied  state  as  a  consequence  of  increased 
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porewater  pressure  and  reduced  effective  stress.  Liquefaction  is 
thus  defined  as  a  changing  of  states  which  is  independent  of 
initiating  disturbance  that  could  be  static,  vibratory,  sea  wave, 
shock  loading,  or  a  change  in  ground  water  pressure.  The 
definition  also  is  independent  of  deformation  or  ground  failure 
movements  that  might  follow  the  transformation,  liquefaction 
always  produces  a  transient  loss  of  shear  resistance  but  does  not 
always  produce  a  long-term  reduction  of  shear  strength**  {Committee 
on  Soil  Dynamics,  1978). 

Liquefaction- “Is  a  phenomenon  where  in  a  mass  of  soil  loses  a 
large  percentage  of  its  shear  resistance,  when  subjected  to 
undrained  monotonic,  cyclic  or  shock  loading,  and  flows  in  a 
manner  resembling  a  liquid  until  the  shear  stresses  acting  on  the 
mass  are  as  low  as  the  reduced  shear  resistance"  (Castro  and 
Poulos,  1977). 

Liquefaction-  “This  term  is  used  to  include  all  phenomena  giving 
rise  to  "a  loss  of  shearing  resistance  or  to  the  development  of 
excessive  strain  as  a  result  of  transient  or  repeated  disturbance 
of  saturated  cohesionless  soils"  (Committee  on  Earthquake 
Engineering,  1985). 

The  definition  by  Seed  (1976)  for  initial  liquefaction  and  that 
for  liquefaction  by  the  Committee  on  Soil  Dynamics  (1978)  will  be  used  for 
this  report.  In  this  way,  blast-induced  liquefaction  can  be  evaluated  by 
comparing  the  amount  of  porewater  pressure  increase  that  occurs  after  the 
passage  of  the  explosive-induced  compressive  stress  wave.  The  porewater 
pressure  increase  above  the  hydrostatic  pressure,  remaining  after  the 
passage  of  the  compressive  stress  wave,  will  be  referred  to  as  the 
residual  excess  porewater  pressure. 


B.  Blast-Induced  Porewater  Pressure  and  Liquefaction 

Documentation  of  blast-induced  liquefaction  exists  in  the 
literature.  Reviews  of  explosive  induced  liquefaction  phenomena  and 
experiences  are  given  by  Blouin  (1978),  Charlie  et  al .  (1980,  1985), 
(Fragaszy  et  al .  (1983)),  Gilbert  (1976),  Marti  (1978),  Rischbieter 
(1977),  Veyera  (1985),  and  others.  Fountains  of  water  and  sand  boils 


-9- 


(Figure  2.1)  have  occurred  and  excess  porewater  pressure  have  been 
measured  following  HE  explosive  tests.  Soils  flowing  towards  and  into 
blast  craters  (Figure  2.2)  have  significantly  changed  the  crater  profile 
and  significant  settlement  outside  the  crater  have  occurred  as  a  result  of 
liquefaction. 

Residual  porewater  pressure  develops  and  liquefaction  occurs 
because  the  stress-strain  behavior  of  granular  soil  is  nonlinear  and 
inelastic.  For  undrained  compressive  loading  of  water  saturated  granular 
soils,  the  soil  skeleton  will  deform  inelastically  while  the  porewater  and 
mineral  phase  (i.e.  soil  grains)  deform  elastically.  Upon  unloading,  some 
or  all  of  the  stress  originally  carried  by  the  soil  skeleton  is 
transferred  to  the  porewater  pressure.  The  development  of  liquefaction 
under  these  circumstances  is  graphically  expressed  in  Figure  2.3  (Fragaszy 
et  al.,  1983). 

A  single  explosion  produces  a  high  frequency,  high  intensity 
compressional  wave  of  very  short  duration.  The  motion  radiates  outward 
from  the  source,  attenuating  with  distance.  The  duration,  magnitude,  and 
form  of  ground  motion  is  a  function  of  the  charge  size  and  shape,  location 
of  the  charge  with  respect  to  the  soil  surface,  and  distance  from  the 
charge  (Lyakhov,  1961).  Compressive  wave  loading  predominates  near  the 
explosive  source  and  shear  and  surface  wave  loading  predominates  at  large 
distances  from  the  explosive  surface  (Melzar,  1978).  For  multiple 
explosive  sources,  the  number  and  timing  of  the  explosive  sources  also 
have  a  large  influence  on  the  duration  and  extent  of  the  induced  ground 
motion. 
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C.  Factors  Influencing  Liquefaction 

Observations  by  many  investigators  both  in  blasting  and  in 

earthquake  research  have  led  to  the  recognition  of  several  significant 

factors  related  to  the  occurrence  of  liquefaction.  These  include: 

Degree  of  Saturation.  The  presence  of  air  In  a  soil  decreases  the 
soil's  potential  for  liquefaction.  Liquefaction  has  been  proven 
difficult  to  induce  in  soils  with  a  degree  of  saturation 
less  than  100  percent  {Florin  and  Ivanov,  1961;  Van  der  Kogel  et 
al . ,  1981;  Studer.,  1977;  True,  1969). 

Gradation.  Particle  Size  and  Shape.  Cohesionless  soils  with  a 
narrow  band  of  gradation  are  generally  considered  easier  to 
liquefy  than  soils  with  a  wide  band  of  gradation  {Damitio,  1978; 
Klohn  et  al.,  1981;  Kummeneje  and  Eide,  1961;  Rischbieter,  1977;). 
The  roundness  of  the  soil  grains  also  seems  to  contribute  to 
liquefaction;  but  angular  sands  have  been  liquefied  under  dynamic 
loading  {Ivanov  et  al.,  1981). 

Compressibility.  The  tendency  for  volume  decrease  is  a  necessary 
condition  for  liquefaction.  Both  the  compressibility  of  the  soil 
skeleton  and  the  fluids  {gas  and  water)  in  the  soil  voids  must  be 
considered  {Florin  and  Ivanov,  1961;  Lyakhov,  1961;  Studer  and 
Kok,  1980;  Van  der  Kogel  et  al.,  1981). 

Permeability  and  Drainage.  The  rate  of  dissipation  of  the 
porewater  pressure  is  a  function  of  the  permeability  of  the  soil, 
the  drainage  length,  and  drainage  boundaries.  Liquefaction  is 
observed  to  occur  mostly  in  soils  which  are  prevented  from 
draining  rapidly  {Florin  and  Ivanov,  1961;  Kok,  1978;  Damitio, 
1978;  Kurzeme,  1971). 

Relative  Density.  Soils  with  a  relative  density  of  less  than  65 
percent  are  typically  considered  susceptible  to  liquefaction 
{Damitio,  1978;  Florin  and  Ivanov,  1961;  Kummeneje  and  Eide,  1961; 
Yamamura  and  Koga,  1974;  Veyera,  1985).  However,  Veyera  (1985) 
and  Hubert  (1986)  have  experimentally  shown  that  shock-induced 
liquefaction  can  also  be  induced  in  soils  at  relative  densities 
above  80  percent. 

Overburden  (Effective  Stress).  The  effective  stress  of  a  soil 
increases  as  the  overburden  increases.  Thus,  the  porewater  must 
reach  a  higher  pressure  to  overcome  the  effective  stress  and  to 
liquefy  the  soil.  Larger  effective  stresses  also  increase  the 
soil  skeleton’s  stiffness,  further  reducing  the  soil's  potential 
for  liquefaction  {Florin  and  Ivanov,  1961;  Rischbieter,  1977; 
Studer  and  Kok,  1980). 
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0.  Laboratory  Studies 

A  summary  of  laboratory  studies  of  shock  induced  liquefaction  is 
given  in  the  following  sections. 

1.  Impact  Loading 

A  simplified  way  of  producing  liquefaction  is  by  impacting  a 
container  filled  with  saturated  sand.  A  primary  drawback  to  this  method 
is  the  lack  of  accuracy  in  estimating  the  impact  force  to  generate 
liquefaction.  Table  2.1  provides  the  summary  of  published  tests  by  impact 
loading. 

2.  Shock  Loading 

Several  approaches  have  been  taken  to  liquefy  a  saturated  material 
by  shock  loading.  Colorado  State  University's  shock  devit-j,  designed  by 
Charlie  and  Veyera  (1985)  under  AFOSR  funding,  uses  water  to  apply  a 
compressive  shock  load.  It  is  the  only  system  in  this  country  that  is 
able  to  vary  the  initial  effective  stress  over  a  wide  range  of  values  and 
record  dynamic  porewater  pressure  responses.  Van  der  Kogel  et  al . ,  1981; 
Pei ry,  1972;  Mason  and  Walter,  1968;  True,  1969;  Studer  and  Hunziker 
(1977)  and  several  others  have  attempted  to  liquefy  saturated  sand  by  air 
blasts.  See  Table  2.2  for  the  summary  of  published  results  of  tests  by 
shock  loading. 

3.  Centrifuge  and  Quasi-Static  Triaxial  Testing 

Schmidt  et  al .  (1981)  and  Fragaszy  et  al .  (1983)  modeled  major 

military  events  by  increasing  the  gravitational  field  in  a  centrifuge. 
The  models  were  subjected  to  scaled  charge  weights  of  the  actual  test  or 
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by  compressed  air.  The  Pacific  Proving  Ground  Tests  on  the  Eniwetok  coral 
atoll  were  reproduced,  resulting  in  liquefaction  (Fragaszy  et  al . ,  1983). 

Fragaszy  and  Voss  (1986)  tested  a  silica  sand  and  the  Eniwetok 
coral  sand  in  a  high  pressure  triaxial  cell.  Several  cycles  of 
quasistatic  isotropic  compressive  loading  produced  liquefaction  at  low 
confining  stresses.  See  Table  2.3  for  the  summary  of  published  results  of 
centrifuge  and  quasi  static  tests. 

E.  Field  Studies 

A  summary  of  field  studies  of  explosive  induced  liquefaction  is 
given  in  the  following  sections. 

1.  Large  Explosives  Tests 

The  detonation  of  up  to  several  kilotons  of  HE  explosions  have 
been  done  by  the  military  to  investigate  ground  motion,  structural 
response  and  crater  formation.  Little  attention  Has  been  given  to  the 
generation  of  porewater  pressure,  except  that  in  some  instances,  the 
observation  of  geysers,  sand  boils,  and  springs  were  observed  and  recorded 
indirectly  indicating  liquefied  layers  at  some  depth  below  the  ground 
water  table  (Langley  et  al . ,  1972;  Banister  et  al . ,  1978;  Banister  and 
Ellett,  1974).  The  location  of  the  ground  water  table  is  ^  primary  factor 
influencing  the  shape  of  the  crater  formed  by  an  explosion.  If  the  scaled 
depth  of  the  ground  water  table  is  very  deep  or  nonexistent,  the  crater 
will  be  in  the  form  of  a  deep  bowl.  In  cohesionless  material  with  a  small 
scaled  distance  to  the  water  table,  the  craters  tend  to  be  very  broad  and 
shallow  which  may  be  caused  by  a  liquefied  state  (Roddy,  1976;  Nordyke, 
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1976;  Melzer,  1978;  Blouin,  1978).  See  Table  2.4  for  the  summary  of 
published  results  for  tests  using  large  HE. 

2.  Explosive- Induced  Compaction 

The  ability  to  compact  loose  saturated  deposits  through  the  use  of 
explosives  has  been  recognized  for  several  decades.  The  stabilization  of 
submerged  slopes  and  construction  sites  by  explosive  densification 
provides  an  economical  method  of  developing  naturally  unsuitable  locations 
in  areas  near  water.  The  general  procedures  for  compacting  large  areas 
involves  simultaneous  detonation  or  delayed  detonation  of  small  buried 
charges  in  a  grid  pattern  (Mitchell  and  Katti,  1981;  Kok,  1978a;  Queiroz 
et  al.,  1967;  Obermeyer,  1980;  Klohn  et  al . ,  1981).  Densification  is 
measured  by  the  amount  of  settlement  in  the  treated  area,  which  can  vary 
from  2  to  10  percent  of  the  total  depth  of  the  cohesionless  soil.  During 
some  of  these  blasting  events,  the  porewater  pressure  response  was 
monitored  (Kummeneje  and  Eide,  1961;  Obermeyer,  1980;  Klohn  et  al . , 
1981).  See  Table  2.5  for  the  summary  of  published  results  of  explosion 
compaction  tests. 

3.  Pits  Away  From  Blast  (Air  Blast  loading) 

Perry  (1972)  and  Rischbieter  et  al .  (1977)  attempted  to  liquefy 
sand  samples  by  placing  them  at  various  distances  from  an  HE  explosive 
event.  The  sand  was  surface  loaded  in  compression  by  the  air  blast.  For 
testing  buoyancy,  Perry  (1972)  placed  objects  of  various  density  buried  at 
different  depths.  Table  2.6  shows  the  summary  of  these  tests. 
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4.  Line  Explosives 

Sanders  (1982)  evaluated  the  liquefaction  potential  as  a  result  of 
a  planned  detonation  by  relating  blast  induced  liquefaction  to  earthquake 
induced  liquefaction.  A  summary  of  the  prediction  is  given  in  Table  2.7. 

5.  Small  Explosives  Tests 

For  charge  weights  of  less  than  1000  kilograms  (most  being  less 
than  10  kilograms),  experiments  have  been  conducted  to  examine  the 
influence  of  various  parameters  on  liquefaction  potential.  The  effects  of 
the  charge  weight,  its  location  relative  to  the  surface,  its  pattern,  and 
its  sequence  of  detonation  (delayed  versus  simultaneous  when  several 
charges  are  involved)  have  been  investigated  by  many  researchers  (Ivanov, 
1967;  Kok,  1977;  Lyakhov,  1961;  Puchkov,  1962;  Rischbieter,  1978;  Dill, 
1966;  Trense,  1977;  Damitio,  1978).  Field  measurements  include  the  peak 
and  residual  porewater  pressure,  the  duration  of  the  excess  porewater 
pressure,  the  settlement  and  occasionally  the  ground  acceleration  at  some 
distance  (Carnes,  1981;  Yamamura  and  Koga,  1974;  Studer  and  Kok,  1980; 
Florin  and  Ivanov,  1961;  Prakash  and  Gupta,  1970;  Arya  et  al . ,  1978; 
Drake,  1978;  Schaepermeier,  1978;  Carnes,  1981).  See  Table  2.8  for  the 
summary  of  published  small  explosive  tests. 

F.  Empirical  Relationships 

1.  Pore  Pressure  Ratio  Based  on  Charge  Weight  and  Distance 

The  porewater  pressure  ratio,  PPR,  also  termed  the  liquefaction 
coefficient,  L,  is  utilized  for  analyzing  the  liquefaction  potential  of  a 
soil.  The  porewater  pressure  ratio  is  the  change  of  the  porewater 
pressure  due  to  loading  related  to  the  initial  effective  stress..  Studer 
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and  Kok  (I960)  developed  the  following  relationship  for  estimating  the 
porewater  pressure  ratio  (liquefaction  coefficient)  from  a  buried 
explosion. 

PPR  =  ur/o^  =  1.65  +  0.64  In  (W1/3  /  R)  Eq.  2.1 

where  PPR  is  the  porewater  pressure  ratio,  ur  is  the  residual  excess 

porewater  pressure,  <j0‘  is  the  initial  effective  stress,  W  is  the  charge 
weight  in  kg,  and  R  is  the  distance  from  the  explosive  in  meters.  The 
porewater  pressure  ratio  (liquefaction  coefficient)  ranges  from  zero  to 
one  with  zero  being  no  increase  in  residual  porewater  pressure  and  one 
being  full  liquefaction.  Equation  2.1  is  plotted  in  Figure  2.4  as  a 
function  of  charge  weight  and  distance.  Figure  2.5  gives  measured  PPR 
from  centrifuge  explosive  tests. 

2.  Liquefaction  Based  on  Charge  Weight 

Lyakhov  (1961)  estimated  that  the  optimum  depth  of  an  explosive,  h 
(in  meters),  for  maximum  radius  liquefaction  is: 

h  =  2.5  W1/3  Eq.  2.2 

where  W  is  the  charge  mass  in  kilograms.  For  an  explosive  charge  buried 
at  the  optimum  depth,  the  maximum  radius  of  liquefaction  can  be  estimated 
from: 

R  =  kx  W1/3  Eq.  2.3 

where  R  is  the  radial  distance  from  the  blast  (in  meters)  and  k  is  an 
empirical  constant  given  in  Table  2.9.  Kok  (1977,  1978)  reported  that  k^ 

has  a  value  of  6.67  (Figure  2.6). 
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Ivanov  (1967)  estimated  that  to  ensure  a  contained  explosion,  the 
charge  weight  per  delay  must  be  less  than: 

W  =  0.055  h3  Eq.  2.4 

where  h  Is  the  burled  charge  depth  In  meters  from  Equation  2.2.  For  a 
fully  contained  explosion,  the  depth  of  liquefaction,  d,  (in  meters)  below 
the  charge  is: 

d  «  1.5  h  Eq.  2.5 

where  h  is  the  buried  charge  depth  In  meters  from  Equation  2.2. 

3.  Pore  Pressure  Ratio  Based  on  Peak  Particle  Velocity,  Strain 
and  Stress 

For  saturated  soils,  Puchkov  (1962)  examined  the  limits  of  the 
peak  particle  velocity  for  liquefaction.  He  concluded  that  a  soil  will 
not  liquefy  when  the  peak  particle  velocity  Is  less  than  7  cm  per  second. 
Obermeyer  (1980)  reported  that  no  significant  increase  in  residual 
porewatev  pressure  were  generated  for  peak  particle  velocities  as  large  as 
2  cm  per  second.  A  safe  maximum  particle  velocity  of  5  to  10  cm  per 
second  is  recommended  for  structural  stability  of  earth  filled  dams 
(Charlie  et  al.,  1985). 

For  single  charges  located  in  water  saturated  soil,  the  peak 
particle  velocity  (in  meters  per  second)  as  a  function  of  the  distance  (in 
meters)  and  the  charge  weight  (in  kilograms)  is  (Drake  and  Ingram,  1981): 

Vpk  =  (7.2)  (R/W1/3)"1,15  Eq.  2.6 

Table  2.10  presents  several  empirical  equations  developed  to 
estimate  peak  stress  in  water  and  water  saturated  soil.  For  a  stress  wave 
traveling  through  an  elastic  media,  a  linear  relationship  exists  on  the 
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wave  front  between  the  peak  compressive  stress,  o^,  and  the  peak  particle 
velocity,  Vp.  This  relationship,  given  by  Kolsky  (1963)  and  Timoshenko 


and  Goodier  (1970)  is: 


=  (P  VJ  V, 


pk 


Eq.  2.7 


where  p  is  the  total  mass  density  of  the  medium  and  Vc  is  the  compressive 


stress  wave  propagation  velocity  through  the  medium.  The  quantity,  "pVc" 

is  commonly  referred  to  as  the  acoustic  impedance  or  specific  acoustic 
resistance  of  the  medium. 

The  peak  compressive  strain,  e  developed  in  a  medium  subjected 

to  a  one-dimensional  compressive  stress  wave  can  be  determined  from  the 
following  equation  (Kolsky,  1963;  Timoshenko  and  Goodier,  1970): 


In  experiments  conducted  by  Veyera  (1985)  on  quartz  sand,  no 
significant  porewater  pressure  increase  was  recorded  for  strains  less  than 
0.005  percent.  Measured  PPR  are  shown  in  Figure  2.7  as  a  function  of  peak 
strain.  The  scatter  in  the  data  points  can  be  explained  by  variations  in 
initial  density  and  initial  effective  stress  of  the  soil.  Liquefaction 
occurred  under  single  compressive  strains  greater  than  0.01  percent  for 
low  effective  stresses  and  low  relative  densities  and  1  percent  for  high 
effective  stresses  and  high  relative  densities.  Utilizing  Equations  2.7 
and  2.8,  for  a  saturated  soil  at  a  void  ratio  of  0.7,  Table  2.11  lists 
relationships  between  peak  compressive  strain,  peak  particle  velocity  and 
a  peak  stress.  A  peak  compressive  strain  of  0.01  percent  corresponds  to  a 
peak  particle  velocity  of  0.15  m/sec  and  a  peak  stress  of  500  kPa.  Veyera 
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(1985)  developed  the  following  relationship  for  estimating  the  porewater 
pressure  ratio  as  a  function  of  peak  strain,  ep)c,  initial  effective 

stress,  0^,  and  Initial  relative  density,  D(,.  The  statistical  best  fit 

of  several  models  examined  is: 

PPR  =  (16.30)  (2epk)*331  (0or,3O8(Drf ,179  Eq.  2.9 

where  Df  and  are  both  in  percent,  and  a ^  is  in  kPa.  Equation  2.9  is 
plotted  in  Figure  2.7  for  Df  equal  to  10  percent,  equal  to  86  k?a,  Dr 
equal  to  80  percent  and  o'Q  equal  to  690  kPa. 

Q.  Analytical  and  Theoretical  Methods 

The  passage  of  stress  waves  associated  with  the  detonation  of  an 
explosive  in  a  saturated  soil  results  in  a  transient  increase  in  porewater 
pressure  and  a  possible  increase  in  residual  excess  porewater  pressure. 
For  saturated  soils,  the  volume  reduction  is  prevented  by  the  increase  in 
porewater  pressure  which,  if  large  enough,  can  produce  liquefaction.  For 
a  period  of  time,  the  liquefied  soil  particles  are  actually  suspended  with 
no  intergranular  stress,  and  the  soil  loses  its  shear  strength.  Following 
liquefaction,  the  soil  consolidates  as  the  particles  settle,  creating  a 
migration  upward  of  excess  porewater.  The  duration  of  the  liquefied  state 
depends  upon  the  drainage  path,  permeability  and  compressibility  of  the 

soil.  As  consolidation  takes  place,  the  soil  regains  its  strength 

(Terzaghi,  1956;  Charlie  et  al . ,  1981;  Blouin  and  Shinn,  1983).  Kim  and 
Blouin  (1984)  used  a  dynamic  finite  element  analysis  to  model  the 

porewater  pressure  during  and  following  the  passage  of  a  compressive 
stress  wave.  Their  analysis  utilized  Biot's  (1962)  theory  of  wave 
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propagation  through  a  porous  media.  Veyera  (1985)  and  Charlie  et  al . 
(1987)  used  a  finite  difference  analysis  also  utilizing  Biot's  theory  to 
analyze  porewater  pressure  Increase  in  laboratory  shock  testing.  These 
analyses  predicted  that  liquefaction  can  be  induced  by  a  compressive 
stress  wave  and  are  in  good  agreement  with  our  experimental  results. 

H.  Summary 

Clear  evidence  of  blast- induced  liquefaction  of  water  saturated 
soils  exist  in  the  literature.  Fountains  of  water  and  sand  boils  have 
occurred  and  excess  porewater  pressure  have  been  measured.  Flow  of  soils 
toward  and  into  craters  and  significant  ground  settlement  surrounding 
craters  have  been  observed.  Empirical  relationships  to  evaluate  blast- 
induced  liquefaction  have  been  developed  from  limited  field  explosive 
tests.  Proper  engineering  judgement  should  be  employed  when  using  such 
relationships  for  specific  conditions.  Only  limited  success  has  been  made 
in  theoretically  understanding  or  analytically  modeling  blast  induced 
porewater  pressure  increases  and  liquefaction.  A  major  problem  in  gaining 
a  better  theoretical  understanding  of  blast  induced  liquefaction  is  the 
complex  stress-strain  behavior  of  granular  material  containing  pore  fluid. 
Soil  is  inherently  a  multiphase  system  and  both  static  and  dynamic  loads 
spoiled  to  a  soil  mass  are  carried  in  part  by  the  mineral  skeleton  and  in 
part  by  the  pore  fluid.  Any  theoretical  or  analytical  model  must 
incorporate  the  response  of  the  mineral  skeleton,  the  pore  fluid  and  their 
interaction. 


Figure  2.1  Liquefaction  is  evident  from  the  geyser ing  of  water  and  sand 
(about  1  m  high)  through  an  instrument  borehole  in  an 
explosion-produced  crater  (Dial  Pack). 


Figure  2.2  Example  of  explosion-induced  liquefaction.  A  large,  300  m 
diameter  crater  was  produced  with  explosives  (Dial  Pack). 


STRESS  (MPa) 


TIME  (min) 


Figure  2.3  Liquefaction  during  the  unloading  phase  from  a  cycle  of 
isotropic  compressive  loading.  Source:  Fragaszy  et  al . ,  (1983) 


CHfiRGE  HEIGHT  (Kg) 


figure  2.4  Liquefaction  coefficient  as  a  function  of  charge  weight  and 
distance  (Studer  and  Kok,  1980). 


(H-U3)/'R  ( (Kg"!/3)/M) 


Figure  2.5  Porewater  pressure  ratio  as  a  function  of  charge  weight  and 
distance  from  explosive  centrifuge  modeling  experiments 
(fragaszy  et  al . ,  1983). 
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Figure  2.6  Liquefaction  coefficent  to  be  expected  as  a  function  of  charge 
and  distance  {Kok,  1977). 
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Figure  2.7  Porewater  pressure  ratio  for  Monterey  No.  Q/3Q  sand  equation 
as  function  of  density  and  effective  stress  (Veyera,  1985). 


Table  2.1  Reference  summary  of  published  research  on  liquefaction  of 
saturated  cohesionless  soil  by  impact  loading. 
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Table  2.2  Reference  summary  of  published  research  on  liquefaction  of 
saturated  cohesionless  soil  by  shock  loading. 
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Table  2.4  Reference  summary  of  published  research  on  liquefaction  of 
saturated  cohesionless  soil  by  large  scale  explosives. 
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Table  2.6  Reference  summary  of  published  research  on  liquefaction  of 
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Table  2.8  continued 
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Table  2.9  Table  of  constant  for  estimated  radius  of  liquefaction  {Ivanov, 
1967). 


Soil* 


Relative  Density 

<%) 


*1 


Fine  Sand 
Fine  Sand 
Fine  Sand 
Medium  Sand 
Medium  Sand 


0-20 
30  -  40 

>  40 
30  -  40 

>  40 


Note  i  Ki  is  for  single  charges 


25  -  15 
9-8 
<_  7 
8-7 
<  6 


Table  2*10  Empirical  Equations  Developed  to  Estimate  Peak  Stress 
in  Water  and  Water  Saturated  Soil 


Reference 

Cole  (1948)* 

°pk 

Lyakhov  (1961) 

°pk 

Drake  and 

Little  (1983) 

°pk 

Crawford  et  al . 

(1974) 

apk 

Equation 

54,90°  (-j^f1*13 
58,900 

Wi/,J 

20,000  {-f^)'2,35 

w 

10,000  (~!^f3-00 


Note:  peak  stress  in  water 

°pk  *  |Jea*(  wa*er  or  porewater  pressure  in  kPa 

R  =  radius  in  meters 
W  =  charge  mass  in  kg 
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Table  2.11  Relationship  Between  Peak  Compressive  Strain, 
Velocity,  and  Peak  Stress  Generated  by  a  Compressive  Shock  Wave  in  a 

Linear  Elastic  Media. 


epk 

..(*>.  - 

^pk 

(m/sec) 

°pk 

(kPa) 

Notes  on  Residual  Porewater 

Pressure 

0.01 

0.15 

500 

no  significant  increase 

0.05 

0.75 

2,500 

liquefaction  of  loose  sand 
under  low  effective  stress 

0.01 

1.50 

5,000 

1.00 

15.00 

50,000 

liquefaction  of  dense  sand 
under  high  effective  stress 

Notes:  saturated  soil  at  void  ratio  =0.7 

low  effective  stress  =  86  kPa 
high  effective  stress  =  690  kPa 
loose  sand  at  10  percent  relative  density 
dense  sand  at  80  percent  relative  density 
epk  s  pea{t  comPressive  strain  in  percent 

Vpk  *  peak  particle  velocity  in  m/sec 

ap(c  =  peak  stress  in  kPa 
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III.  LABORATORY  SHOCK  TESTS  ON  SAND 

The  results  of  laboratory  shock  tests  which  we  conducted  on  five 
saturated  cohesionless  soils  are  presented  in  this  chapter.  An  overview 
of  our  equipment,  instrumentation,  test  procedures  is  also  presented. 

A.  Introduction 

The  one-dimensional,  shock  loading  system,  designed  by  Charlie  and 
Veyera  (1985),  was  used  to  investigate  the  porewater  pressure  response  of 
five  saturated  sands  subjected  to  shock  loading.  The  shock  apparatus, 
funded  by  the  United  States  Air  Force  Office  of  Scientific  Research  (Grant 
AFOSR-80-026G) ,  is  located  in  the  Geotechnical  Engineering  Research 
Laboratory  of  the  Civil  Engineering  Department  of  Colorado  State 
University.  For  greater  details  on  the  equipment,  instrumentation  and 
procedures,  see  Veyera  (1985)  and  Charlie  et  al .  (1985). 

B.  Test  Equipment  and  Instrumentation 

To  simulate  an  explosive  loading  condition,  high  amplitude 
compressive  stress  waves  with  submi 11 i second  rise  time  to  peak  are  needed. 
To  meet  this  requirement,  the  shock  loading  system  is  designed  to  apply  a 
compressive  shock  wave  loading  of  up  to  35,000  KPa  to  a  saturated  sand 
sample.  Cross  sectional  views  of  the  sample  container,  confining  pressure 
tube,  piston,  momentum  trap  and  the  location  of  the  porewater  pressure 
transducers  and  pressure  ports  and  valves  are  shown  in  Figure  3.1. 


Transient  and  static  porewater  pressure  was  measured  with 
porewater  pressure  transducers.  The  piezoresisti ve,  silicon  diaphragm, 
strain  gauge  pressure  transducers  are  air  blast  transducers  modified  by 
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the  manufacturer  by  mounting  a  small  perforated  plate  over  the 
transducer's  casing  (Endevco  Model  8511a-5KMl).  This  modification  allowed 
the  static  and  transient  porewater  pressure  to  be  measured  in  the  sample. 

A  dual  time  base  digital  waveform  recorder  was  used  for  recording 
the  porewater  pressure  at  two  sampling  rates:  rapid  at  the  beginning  as 
the  stress  wave  passes  to  record  the  transient  porewater  pressure  response 
and  then  slower  to  record  the  residual  pore  pressure.  The  computer  system 
to  store,  manipulate  and  analyze  the  data  consisted  of  a  desktop  computer, 
a  floppy  disk  device  for  mass  storage,  a  dot-matrix  line  printer,  and  a 
graphic  plotter  (Figure  3.2). 

Figure  3.3  shows  the  shock  facility  prepared  for  loading.  The 
photograph  shows  the  monentum  trap  in  the  foreground,  followed  by  the 
sample  container,  the  confining  pressure  tube  and  the  air  cannon  the 
background. 

C.  Physical  Properties  of  the  Soils 

Several  index  tests  were  performed  on  each  sample  to  investigate 
various  physical  properties.  All  tests  were  conducted  according  to  the 
standard  laboratory  procedures  where  applicable  and  include  the  following: 
grain  size  analysis  (ASTM  0422), 

-  soil  classification  (ASTM  D2487), 
specific  gravity  (ASTM  0854), 

-  relative  density  determination  (ASTM  D2049), 
photomicrograph  (Bureau  of  Reclamation), 


spectrographic  analysis  (Bureau  of  Reclamation),  and 
skeleton  stress-strain  curves  (Hendron,  1963) 


Test  results  are  given  in  later  sections  of  this  report. 

The  water  saturated  sand  samples  used  in  this  experimental 
investigation  represent  a  two-phase  medium.  The  presence  of  the  solid 
particles  in  the  water  affects  the  density  and  the  compressibility  of  the 
soil-water  mixture.  These  changes  in  density  and  compressibility  need  to 
be  considered  in  evaluating  compressive  stress  wave  propagation  velocity, 
peak  stress  and  peak  strain.  A  procedure  presented  by  Ri chart  et  al . 
(1970)  was  followed  and  will  be  outlined  here. 

The  total  mass  density,  pt,  of  a  fully  saturated  soil -water 

mixture  can  be  determined  from: 

G  +  e 

pt  =  pw"n^r  Eq- 

where  G  is  the  specific  gravity  of  the  solid  particles,  p  is  the  mass 
s  w 

density  of  water  and  e  is  the  void  ratio  defined  as  the  volume  of  voids 

divided  by  the  volume  of  solids. 

The  compressibility  of  a  fully  saturated  soil-water  mixture, 
considering  the  solid  particles  to  be  suspended  in  the  water,  consists  of 
two  contributing  factors:  the  compressibility  of  the  solid  particles  and 
the  compressibility  of  the  fluid.  The  mixture  compressibility  can  be 
determined  from: 

-r1— =  r~  i-rfr>  *  r~  <-rrr>  ^  3-2 

mix  w  s 
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where  B  Is  the  bulk  modulus  of  the  water  and  B  is  the  bulk  modulus  of 
w  s 

the  soil  particles.  For  quartz  particles,  Is  about  30,680  MPa  and  for 

distilled,  de-aired  fresh  water  at  20  degrees  Celsius,  B  is  about  2,140 

w 

MPa  {Richart  et  al ,  1970).  It  should  be  noted  that  the  theory  of  mixtures 
assumes  that  the  solid  particles  are  suspended  in  the  water.  The  total 
stress  will  be  somewhat  larger  since  the  solid  particles  are  actually  in 
contact  with  one  another.  For  small  strain  conditions,  this  difference  is 
small  and  can  neglected. 

Using  the  bulk  modulus  and  the  total  mass  density  for  the  mixture, 

the  compressive  stress  wave  propagation  velocity  through  the  mixture, 

V  .  ,  can  be  found  from: 
mix 


B  . 

V  -  r  mix  -iG.5 
mix  “  L  pt  j 


Eq.  3.3 


The  value  of  Vm-x  includes  the  effects  of  density  and  compressibility. 

Calculated  values  of  V  •  for  quartz  sand  (G  =  2.65)  and  coral  sand  (G  = 

im  x  s  s 

2.80)  are  given  in  Table  3.1.  The  values  for  V  .  for  the  granitic  sand 

ill  1  A 

and  gravels  are  very  similar  to  those  for  quartz  sand.  By  substituting 
V  .  for  V  and  p.  for  p,  Equations  2.7  and  2.8  can  also  be  used  to  relate 

III  I  A  v  w 

peak  particle  velocity,  peak  strain  and  peak  stress  in  ?  soil-water 
mixture. 


0.  Variation  of  Parameters 

The  following  six  parameters  were  varied  in  a  systematic  and 
controlled  manner  for  this  study: 
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-  the  initial  relative  density  of  the  samples, 

-  the  initial  effective  stress  on  the  samples, 

-  the  shape  of  the  soils  grains  {subrounded  and  angular) 

-  the  size  of  the  soil  grains  (silt,  sand  and  gravel) 

the  hardness  of  the  soil  grains  {quartz,  granite  and 
coral ) 

the  number  and  intensity  of  applied  shock  loads. 

The  Initial  relative  density  of  the  samples  used  in  the 
experimental  investigations  ranged  from  0  to  100  percent.  To  obtain  a 
sample  near  zero  relative  density,  the  sand  was  placed  using  a  metal 
funnel.  The  funnel  was  slowly  lifted  upward  and  simultaneously  rotated 
inside  the  sample  container  while  keeping  the  spout  about  1.27  cm  above 
the  placed  material.  The  preparation  of  samples  at  relative  densities 
greater  than  zero  percent  and  less  than  one  hundred  percent  was 
accomplished  using  a  method  presented  by  Ladd  (1978)  called  the 
undercompaction  method.  Ten  layers  of  equal  weight  were  placed  using  the 
funnel  approach  described  previously.  The  layers  were  individually 
compacted  to  successively  higher  percentages  of  the  final  sample  density, 
varying  linearly  by  layer.  The  first  layer  was  five  percent 
undercompacted.  The  preparation  of  samples  at  relative  densities  of  one 
hundred  percent  was  accomplished  by  funnel  placement  followed  by  vibration 
of  the  sample  container.  Soils  were  subjected  tc  initial  effective 
stresses  ranging  from  86  kPa  to  690  kPa.  The  sample's  back  pressure  was 
maintained  at  3A5  kPa  with  the  confining  pressure  adjusted  according  to 
the  required  effective  stress.  The  samples  were  loaded  by  compressive 
stress  wave  loadings  between  0.10  MPa  and  4.2  MPa  per  impact. 
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To  evaluate  the  influence  of  each  factor,  a  statistical  analysis 
of  the  data  was  performed  to  develop  3  model  for  predicting  the  porewater 
pressure  ratio,  PPR  (also  termed  the  liquefaction  coefficient,  L).  The 
porewater  pressure  ratio  is  a  nondimensionalized  factor  defined  as: 

u„ 


PPR  = 


Eq.  3.4 


where 

cr‘  =  the  initial  effective  stress, 
o 

ur  =  the  residual  change  in  porewater  pressure  after  passage  of 
the  stress  wave. 

In  nondimensionalized  form,  the  changes  in  porewater  pressure  following 
loading  of  soils  having  differing  initial  effective  stress  can  be 
compared.  The  results  from  laboratory  and  field  investigations  can  also 
be  compared. 


E-  Sample  Saturation 

To  ensure  full  saturation  of  the  sample,  carbon  dioxide  gas  was 
introduced  at  the  sample's  bottom  and  then  distiled,  de-aired  water  was 
slowly  introduced.  After  about  three  pore  v  ;umes  of  water  had  passed 
through  the  sample,  a  backpressure  Of  345  kPa  was  applied  and  the  sample's 
porewater  pressure  response  was  "quasistatically"  checked  to  date* ..line  the 
degree  of  saturation.  It  was  observed  that  the  porewater  pressure 
response  decreased  with  increasing  sample  density  and  effective  stress. 
Calculations  of  the  porewater  pressure  ratio  for  saturated  sands  utilizing 
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methods  suggested  by  Skempton  (1954)  using  skeleton  stress-strain 

relationships,  showed  similar  trends.  The  C-parameter  is  defined  as 

C  =  Au/Ao  Eq.  3.5 

where  au  is  the  change  in  porewater  pressure  and  Ao  is  the  change  in 

confining  pressure.  Lambe  and  Whitman  (1969)  present  the  following 

relationship  to  calculate  C,  where  0  is  the  constrained  modulus  of  the 

soil  skeleton  and  B  .  is  the  bulk  modulus  of  the  soil  particle  and  water 

mix 

mixture,  and  is  given  by  Equation  3.2. 

C  =  1/(1  +  D/Bm-X)  Eq.  3.6 

Utilizing  Equations  3.2  and  3.6  and  laboratory  tests  to  obtain  the 
constrained  modulus,  the  theoretical  C-parameter  for  100  percent 
saturation  may  be  found.  Increases  in  the  contrained  modulus  of  the  soil 
skeleton  will  affect  the  C-parameter  by  lowering  its  value  below  1.0  (Lee 
et  al.,  1969).  When  the  C-parameter  was  close  to  its  theoretical  value, 
the  soil  was  considered  saturated.  Measured  compressive  stress  wave 
propagation  velocities  through  the  samples  also  showed  they  were  saturated 
at  C  values  close  to  that  given  by  Equation  3.6. 


Shock  testing  was  conducted  by  Veyera  (1985)  to  evaluate  residual 
porewater  pressure  Increases  and  liquefaction  potential  of  Monterey  No. 
0/30  sand  as  a  function  of  peak  shock  induced  compressive  strain,  soil 
density,  and  effective  stress.  Details  of  the  grain  size  analysis, 
classification  and  relative  density  tests  are  given  by  Muzzy  (1983)  and 
Charlie  et  al.  (1985).  The  soil  was  tested  at  four  relative  densities  of 
0,  20,  40  and  80  percent  under  four  effective  stresses  of  86,  172,  345  and 
690  kPa.  The  sand's  characteristics  are  summarized  below  and  in  Appendix 
A.  2. 

1.  Grain  Size  Distribution  and  Grain  Shape 

Monterey  No.  0/30  is  a  fine,  uniform,  poorly  graded  sand  ($P)  with 
less  than  one  percent  of  the  material  being  finer  than  0.150  mm.  The  Dgg 

particle  size  is  0.45  mm.  The  coefficients  of  uniformity  and  curvature 
are  1.65  and  1.00,  respectively.  The  grain  shape  is  subrounded  to 
subangular. 

2.  Shock  Induced  Porewater  Pressure  Ratio 

To  evaluate  the  influence  of  initial  effective  stress,  initial 
sample  density  and  applied  compressive  strain,  a  statistical  multivariate 
regression  analysis  of  the  data  was  performed  to  develop  a  model  for 
predicting  the  porewater  pressure  ratio,  PPR,  defined  in  Equation  3.4.  A 
summary  of  numerical  results  of  all  samples  tested  is  given  in  Tables  A. 2 
through  A. 12. 
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The  statistical  best  fit  of  all  the  models  examined  is: 

Eq.  3.7 


PPR  =  (16.30)  (£cpk)*331  (o;)**308(Dr)'-179 


where  Dr  and  epk  are  both  in  percent,  and  is  in  kPa.  The  £epk  quantity 

is  the  cumulative  compressive  str. in  of  the  sample.  For  single  (first) 
impacts,  it  represents  the  peak  strain.  For  multiple  strains  under 
undrained  conditions,  the  term  represents  the  peak  strain  applied  at 

a  given  loading  plus  the  peak  strain  induced  during  previous  impacts.  The 

coefficient  of  determination  (R  )  and  the  standard  error  of  estimate  (S) 
for  Equation  3.7  are  65.9  percent  and  0.187,  respectively. 

The  statistical  model  considering  only  the  data  from  the  first 
impact  loadings  from  each  data  set  has  the  following  form: 

PPR  =  (16.00)  Upk)°-331  (a;f°-308  (DRf°-179  Eq.  3.8 

Equation  3.7  is  plotted  with  the  experimental  PPR  vs  peak  strain  in  Figure 


3.4. 
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S.  Quartz  Beach  Sand  (Tyndall  AFB,  Florida) 

The  Tyndall  beach  sand  was  obtained  at  a  location  500  m  east  of 
the  NATO  blast  facility  at  Tyndall  AF8,  Panama  City,  Florida.  The  Tyndall 
beach  sand,  a  fine,  subrounded,  quartz  sand,  was  shock  tested  by  Hubert 
(1986)  at  relative  densities  of  55,  63  and  at  73  percent  under  an 
effective  stress  of  345  kPa.  The  sand's  characteristics  are  given  below 
and  summarized  in  Table  A. 14  in  Appendix  A. 2. 


1.  Grain  Size  Distribution  and  Grain  Shape 

The  soil  is  a  uniform,  poorly  graded  material  (SP)  with  a  size 

of  0.22  mm.  The  minimum  particle  size  is  greater  than  0.075  mm  (#200 
sieve).  The  grain  shape  is  subrounded.  The  coefficients  of  uniformity 
and  curvature  are  1.47  and  1.04,  respectively. 


2.  Shock  Induced  Porewater  Pressure  Ratio 

A  multivariate,  regression  analysis  was  not  performed  due  to  a 
limited  amount  of  sand  (only  3  tests  could  be  performed).  Figure  3.5 
shows  Veyera's  (1985)  equation  for  Monterey  No.  0/30  sand  relative  to  the 
data  for  the  55,  63  and  75  percent  relative  density  tests  under  an 
effective  stress  of  345  kPa.  Veyera's  (1985)  equation  for  medium  quartz 
beach  sand  fits  the  actual  data  for  the  Florida  sand  quite  well.  Table 
A.  15  presents  the  results  of  the  shock  tests. 


-45- 


H.  Granitic  River  Sand  (Poudre  River,  Colorado} 

The  Poudre  Valley  Sand,  an  angular  granitic  sand,  was  shock  tested 
by  Hubert  (1986).  Nine  samples  of  the  Poudre  Valley  sand  were  impacted 
until  liquefaction  occurred  (porewater  pressure  ratio  reached  1.0),  The 
soil  was  tested  at  four  relative  densities  of  20,  40,  60  and  80  percent 
under  three  effective  stresses  of  86,  172,  and  517  kPa.  The  sand 
characteristics  are  summarized  below  and  given  in  Table  A. 16  in  Appendix 
A. 3. 

1.  Grain  Size  Distribution  and  Grain  Size 

The  sand  is  poorly  graded  (SP)  with  a  DgQ  size  of  0.52  mm.  The 

grain  shape  is  angular.  The  coefficients  of  uniformity  and  curvature  are 
3.85  and  0.95,  respectively. 

2.  Shock  Induced  Porewater  Pressure  Ratio 

The  results  were  evaluated  by  using  linear  regression  analysis. 
The  best  fitted  equation  for  the  data  is: 

PPR  =  10.59  (Sepk)'429  (o)0''171  (Drf '181  Eq.  3.9 

The  coefficient  of  determination  (R2)  and  the  standard  error  of  estimate 
(S)  are  61.0%  and  0.2114  respectively. 

In  Table  A. 17,  the  results  of  the  nine  tests  may  be  found.  In 
Figure  3.6,  Equation  3.9  is  plotted  with  the  combined  data. 

Most  tests  were  found  to  liquefy  upon  the  first  impact  when  peak 
compressive  strains  exceeded  0.2  percent.  Liquefaction  could  not  be 
induced  even  under  multiple  impacts  if  the  strains  were  less  than  0.01 
percent.  Between  0.01  and  0.2  percent  strain,  the  number  of  impacts 
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required  to  induce  liquefaction  was  found  to  be  a  factor  of  strain 
amplitude,  initial  effective  stress  and  initial  relative  density. 

I.  Granitic  River  Fine  Sand  and  Gravel  (Poudre  Valley,  Colorado) 

Shock  testing  was  conducted  by  Chouicha  (1987)  to  evaluate  the 
effects  of  grain  size  (fine  sand  vs  gravel)  on  the  liquefaction  potential. 
Tests  were  conducted  on  samples  of  Poudre  Valley  sand  and  gravel  sieved 
and  recombined  to  give  the  desired  grain  size  distribution.  Each  soil  was 
tested  at  relative  densities  of  30  and  70  percent  under  effective  stresses 
of  207  and  345  kPa.  The  fine  sand  and  gravel  characteristics  are 
summarized  below  and  given  In  Table  A. 18  in  Appendix  A. 4. 

1.  Grain  Size  Distribution,  Fine  Sand 

The  sand  is  uniform,  poorly  graded  (SP)  with  a  minimum  particle 
size  greater  than  0.075  mm  (#200  sieve).  The  grain  shape  is  angular,  the 
050  particle  size  for  this  material  is  0.20  mm.  The  coefficients  of 

uniformity  and  curvature  are  1.50  and  1.11,  respectively. 

2.  Grain  Size  Distribution,  Gravel 

The  gravel  is  uniform,  poorly  graded  soil  (GP)  with  less  than  4 
percent  of  the  material  being  finer  than  4.75  mm  (#4)  and  the  D^g  particle 

size  is  8.9  mm.  The  grain  shape  is  angular.  The  coefficients  of 
uniformity  and  curvature  are  1,50  and  1.07,  respectively. 
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3.  Shock  Induced  Porewater  Pressure  Ratio 

Figures  3.7  and  3.8  show  the  effect  of  mean  grain  size  on  the 
porewater  pressure  ratio  and  the  liquefaction  potential  for  the  fine  sand, 
the  gravel  and  for  data  from  both  combined.  The  general  trend  is  for 
gravel  to  be  more  resistant  to  liquefaction.  In  other  words,  the  pore 
pressure  ratio  increases  with  decreasing  mean  grain  size  at  a  given  sample 
sum  of  the  peak  strains.  The  test  results  are  given  in  Tables  A. 19  to 
A. 22. 

The  best  fit  statistical  models  have  the  following  forms: 

Fine  Sand  Data 


PPR  «  114.82  (Eepk)0-518  (o;f0-329  <Drf°-380  Eq.  3.10 
(R2  =  94.3  percent  and  S  *  0.0353) 

Gravel  Data 


PPR  =  15.85  (2cpk)0'400  (a;)*0'236  (Dr)**°* 119  Eq.  3.11 
(R2  =  92.9  percent  and  S  =  0.0329) 


All  Data  (fine  sand  and  gravel) 


PPR  =  34.67  (^pk)°'44  (o;f°-264  (Drf0,236  Eq.  3.12 


(R2  =  91.2  percent  and  S  =  0,0387) 
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Equations  3.10,  3.11  and  3.12  are  plotted  on  Figures  3.7,  3.8  and  3.9, 
respectively. 

J.  Coral  Beach  Sand  (Enlwetok,  South  Pacific) 

Shock  testing  was  conducted  by  Hubert  (1986)  to  evaluate  the 
residual  porewater  pressure  increases  and  liquefaction  potential  of  coral 
sand.  The  soil  was  tested  under  relative  densities  of  0,  20,  40,  60,  80 
and  100  percent  under  effective  stresses  of  172,  345,  and  517  kPa.  The 
Eniwetok  coral  sand  was  furnished  by  the  Air  Force  Weapons  Laboratory  in 
Albuquerque,  New  Mexico  and  by  the  U.S.  Army  Waterways  Experiment  Station 
in  Vicksburg,  Mississippi.  The  coral  sand  characteristics  are  summarized 
below  and  given  in  Table  A. 23  in  Appendix  A. 5. 

1.  Grain  Size  Distribution  and  Grain  Shape 

The  physical  properties  of  the  sand  were  included  with  the  sand 
shipment  and  were  used  for  this  experiment.  The  grain  size  analysis 
indicates  a  uniformly  graded  sand  (SP),  Dgg  equal  to  0.48  mm  and  no 

particles  smaller  than  the  number  200  sieve  (0.075  mm).  The  grain  shape 
is  subrounded  to  subangular.  The  coefficients  of  uniformity  and  curvature 
are  1.66  and  1.09  respectively. 


2.  Induced  Porewater  Pressure  Ratio 

A  summary  of  the  tests  results  are  given  in  Tables  A, 24  to  A. 29 
and  plotted  in  Figure  3.10.  The  best  fit  data  for  the  PPR  is: 


PPR  =  5.81  {Zc 


pk' 


,0.429  ,  -0.176  ,n  ,-0.022 

'  (n  1  (Dd) 


<°0> 


with 


equal  to  61.6  percent  and  S  equal  to  0.272. 


Eq.  3.13 
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|  K.  Summary 

[r  ■  Four  saturated  cohesionless  sands  and  one  gravel  were  subjected  to 

|  compressive  shock  wave  loading.  The  best  fit  multivariate  regression 

| 

model  to  predict  the  PPR  as  a  function  of  peak  strain,  effective  stress 
i  and  relative  density  for  each  soil  tested  are: 

y 

!  Monterey  No.  0/30  (Subrounded  Quartz  Sand,  Dgg  =  0.45  mm) 

PPR  *  (16.30)  Uepk)*331  (a;)"’308  (Dr)--179  Eq.  3.7 

Tyndall  AFB  Sand  (Subrounded  Quartz  Sand,  DgQ  -  0.22  mm) 

(multivariate  regression  not  conducted  -  Eq.  3.7  fits  data) 

Poudre  Valley  Sand  (Angular  Granitic  Sand,  =  0.52  mm) 

PPR  =  (10.59)  (£epk)'429  (o^)"’171  (°R)"'181  Eq.  3-9 

Poudre  Valley  Fine  Sand  (Angular  Granitic  Sand,  D^g  =  0.20  mm) 


PPR  =  (114.82)  Uepk)'518  (a^f*329  (DR)''380  Eq.  3.10 


3 
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Poudre  Valley  Gravel  (Subangular  Granitic  Sand,  DgQ  ~  0.90  mm) 


PPR  =  (15.85)  <zepk)‘400  (a;)-236  (DRf  *119  Eq.  3.11 


Coral  Sand  (Subrounded  Coral  Sand,  D^g  =  0.48  mm) 


PPR  =  (5.81)  Uepkr429  (a‘)~A76  (DRf -022 


Eq.  3.13 


where  epk  ancl  ar8  *n  Percent  aric*  ao  ™  kpa* 

To  compare  these  multivariate  regression  models,  the  PPR  predicted 
by  each  of  these  equations  are  plotted  on  Figure  3.11  as  a  function  of 
peak  strain  for  0R  =  50%  and  a °  =  100  kPa.  All  equations  predict 

liquefaction  (PPR  =  1)  when  the  peak  strain,  epk,  exceeds  0.1  percent. 
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.  1 ( a  and  b)  Cross-section  of  the  confining  pressure  tube  and  the 
sample  container  (Veyera,  1985). 


Figure  3.3  Experimental  shock  facility  prepared  for  loading  (view  from 
momentum  trap). 
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Figure  3.4  Porewater  pressure  ratio  for  Monterey  No. 

function  of  density  arid  effective  stress. 
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figure  3.5  Porewater  pressure  ratio  for  Tyndall  beach  sand  as  a  function 
of  the  strain. 
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Figure  3.9 


Porewater  pressure  ratio  as  a  function  of  the  sum  of  the 
strains,  Poudre  Valley  fine  sand  and  gravel. 


peak 


Figure  3.10  Porewater  pressure  ratio  for  Eniwetok  coral  sand  as  a 
function  of  the  strain,  relative  density  and  effective 
stress. 
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Honterey  No.  0/30  Sand  Eq  3.7 

Poudre  Valley  Sand  Eq  3-9 

Poudre  Valley  Fine  Sand  Eq  3.10 
Poudre  Valley  Gravel  Eq  3-11 

Eniwetok  Coral  Sand  Eq  3.13 
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Figure  3.11  Porewater  pressure  ratio  as  a  function  of  the  strain  for  the 
five  soils  tested  (0R  =  50%  and  o'  =  100  kPa). 


Table  3.1  (a)  Stress  Wave  Propagation  Parameters  for 
Quartz  Sand  (6S  =  2.65) 


« 

(Kg/M3) 

B»i* 

(KPa) 

V.ix 

V  .  /V 
nix  « 

Vc 

(Kg/(M2-seo) 

.803 

1915 

4421407 

1519 

1.013 

2908885 

.779 

1926 

4485912 

1526 

1.017 

2940602 

.755 

1940 

4554171 

1532 

1.021 

2972080 

.731 

1953 

4626521 

1539 

1.026 

3005667 

.707 

1967 

4703343 

1546 

1.031 

3040982 

.683 

1980 

47  85064 

1555 

1.037 

3078900 

.659 

1995 

4872167 

1563 

1.042 

3118185 

.635 

2009 

4965202 

1572 

1.048 

3158148 

.611 

2024 

5064797 

1582 

1.055 

3201968 

.587 

2040 

5171672 

1592 

1.061 

3247680 

.563 

2056 

5286654 

1604 

1,069 

3297824 

Table  3.1  (b)  Stress  Wave  Propagation  Parameters  for 
Coral  Sand  (G$  3  2.80) 


mix 


mix 


V  .  /V 
mix  w 


(Mg/m3) 

(MPa) 

(m/s) 

.818 

1.541 

4579 

1724 

1.149 

.797 

1.560 

4640 

1725 

1.150 

.777 

1.579 

4702 

1726 

1.151 

.755 

1.598 

4774 

1728 

1.152 

.735 

1.617 

4843 

1730 

1.153 

.713 

1.637 

4922 

1734 

1.156 

.693 

1.657 

4999 

1737 

1.158 

.673 

1.676 

5080 

1741 

1.161 

.652 

1.698 

5170 

1745 

1.163 

.631 

1.720 

5267 

1750 

1.167 

IV.  LABORATORY  SHOCK  AND  QUASI-STATIC  TESTS  ON  SILT 


This  chapter  summarizes  our  results  of  uniaxial  quasi-static  and 
shock  loading  laboratory  tests  conducted  by  Bolton  (1988)  on  saturated 
quartz  beach  sand,  on  saturated  clayey  silt  and  on  saturated  sand  clayey 
silt  mixtures.  These  tests  were  run  at  the  request  of  DNA. 

A.  Introduction 

Very  little  information  is  available  about  the  porewater  pressure 
response  of  water  saturated  silts  subjected  to  explosive  induced  shock 
waves.  To  gain  needed  data,  a  series  of  uniaxial  laboratory  quasi -static 
and  shock  loading  tests  were  conducted  on  a  sand,  on  a  clayey  silt,  and  on 
a  sand-clayey  silt  mixture.  Tests  were  conducted  at  various  initial  void 
ratios  and  initial  effective  stresses.  The  quasi-static  tests  subjected 
the  samples  to  slowly  increasing  confining  stress  (no  inertia  effects)  and 
the  resultant  increase  in  the  sample's  porewater  pressure  was  recorded. 
Following  loading  the  confining  stress  was  slowly  reduced  to  the  original 
value  and  the  decrease  in  porewater  pressure  was  recorded.  The  excess 
porewater  pressure  after  unloading  is  termed  the  quasi -static  residual 
porewater  pressure.  The  shock  loading  tests  subjected  the  sample  to  a 
transient  compressive  shock  wave  having  submill isecond  rise  time  to  peak 
and  the  resultant  increase  in  the  sample's  porewater  pressure  was 
recorded.  The  excess  porewater  pressure  after  the  stress  wave  has  passed 
is  termed  the  shock  induced  residual  porewater  pressure. 


8.  Test  Equipment  and  Procedure 

The  one-dimensional,  confined  loading  system,  designed  by  Charlie, 
Veyera  and  Huzzy  (1982)  and  used  by  Charlie,  et  al .  (1985),  was  utilized 
by  Bolton  (1988)  for  this  research.  Details  of  the  experimental  equipment 
and  instrumentation,  as  well  as  the  experimental  set-up  and  testing 
procedures,  are  discussed  briefly  in  Chapter  III  of  this  report.  Full 
details  can  be  found  in  Charlie  et  al.  (1985). 

C.  Test  Instrumentation 

To  measure  the  porewater  pressure  time  history,  special  oorewater 
pressure  transducers  and  high  speed  recording  systems  were  utilized. 
These  components  are  discussed  briefly  in  the  following  sections. 

1.  The  Porewater  Pressure  Transducers 

Pi ezoresi stive  porewater  pressure  transducers  were  used  to  measure 
quasi-static  transient  and  static  porewater  pressure.  The  piezoresistive, 
silicon  diaphragm,  strain  gauge  pressure  transducers  are  air  blast 
transducers  modified  by  the  manufacturer  by  mounting  a  stainless  steel 
perforated  plate  over  the  transducer's  pressure  sensor  (ENDEVCO  Model 
8511a-5KMl)  .  This  modification  allowed  the  sample's  porewater  pressure 
to  be  measured. 

2.  Data  Recording  System 

For  the  quasi-static  tests,  a  strip-chart  recorder  was  utilized  to 
record  both  the  applied  pressure  and  the  sample's  porewater  pressure 
response.  For  the  shock  tests,  a  dual  time  base  digital  waveform  recorder 
was  used  for  recording  the  applied  pressure  the  the  sample's  porewater 
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pressure  at  two  sampling  rates,  rapid  at  the  beginning  as  the  stress  wave 
passes  to  record  the  transient  pcrewater  pressure  response  and  then  slow 
to  record  the  residual  pore  pressure  were  used.  The  computer  system  to 
store,  manipulate  and  analyze  the  data  consisted  of  a  desktop  computer,  a 
floppy  disk  device  for  mass  storage,  a  dot-matrix  line  printer,  and  a 
graphics  plotter. 

0.  Physical  and  Index  Properties  of  the  Soils 

Basic  physical  and  index  properties  for  Monterey  No.  0/30  sand. 
Bonny  silt  and  the  sand-silt  mixture  are  given  in  Table  4.1.  The  sand- 
silt  mixture  consisted  of  equal  weights  of  Monterey  No.  0/30  sand  and 
Bonny  silt.  Monterey  No.  0/30  sand  is  a  poorly  graded,  sub-rounded, 
quartz  beach  sand  (SP)  obtained  from  Monterey,  California  and  Bonny  silt 
is  a  wind  blown  deposit  of  clayey  silt  (MC)  and  is  typical  of  loess 
deposits  found  in  Colorado,  Kansas,  Nebraska  and  Wyoming.  The  Bonny  silt 
was  obtained  by  the  Bureau  of  Reclamation  from  the  borrow  area  at  Bonny 
Reservoir  in  eastern  Colorado. 

E,  Quasi-Static  Tests 

Samples  of  sand  and  sand-clayey  silt  mixtures  were  placed  at  void 
ratios  of  0.66  and  0.80  and  samples  of  clayey  silt  were  placed  at  a  void 
ratio  of  0.95.  Figure  4.1  presents  the  gradation  for  these  three  soils. 
These  samples  were  then  consolidated  to  an  initial  effective  confining 
stress  of  172  kPa.  The  air  dried  soil  was  placed  and  compacted  in  ten 
layers.  To  ensure  uniform  soil  density,  the  under  compaction  method  was 
utilized  with  the  first  layer  being  5  percent  undercompacted  and  the  last 
layer  compacted  to  the  required  void  ratio.  To  ensure  full  saturation. 
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the  soil  was  first  flushed  with  carbon  dioxide  gas,  a  vacuum  was  applied, 
then  three  pore  volumes  of  deaired  water  was  introduced  info  the  bottom  ot 
the  sample,  and  a  backpressure  of  345  kPa  was  applied  to  the  porewater. 
Table  4.2  lists  the  symbols  used  for  the  tests  listed  in  the  following 
sections. 

1.  Sand 

The  residual  porewater  pressure  results  from  the  quasi -static 
tests  on  the  saturated  sand  at  void  ratios  of  0.66  and  0.80  are  given  in 
Tables  4.3  and  4.4  Figures  4.3  and  4.4.  Figures  4.5,  4.6  and  4.7  present 
the  residual  porewater  pressure  increase  and  PPR  as  a  function  of  peak 
strain.  The  best  fit  multivariate  regression  models  are: 

PPR  =  60  (Sepk)0'41  (a;f-51  (DR)*14  Eq.  4.1 
PPR  =  43  (£epk)°‘41  (°;)~'51  Ed-  4.2 

where  epk  is  percent  strain,  op  is  the  effective  stress  in  kPa,  e  is  the 
void  ratio  and  Dr  is  the  relative  density  in  percent. 

2.  Sand-Clayey  Silt  Mixture 

The  results  of  the  quasi -static  tests  on  the  saturated  Monterey 
No.  0/30  sand  and  8onny  silt  mixture  at  void  ratios  of  0.66  are  given  in 

Table  4.5  and  Figures  4.8  and  4.9.  The  best  fit  mul ti-regression  model  is 

given  in  Equation  4.3.  The  sand-clayey  silt  mixture  and  the  clayey  silt 

were  analyzed  together  since  there  were  no  dramatic  difference  in 


behavior. 
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3.  Clayey  Silt 

The  results  of  the  quasi-static  tests  on  Bonny  silt  at  a  void 
ratio  of  0.95  are  given  in  Table  4.6  and  Figures  4.8  and  4.9.  The  best 
fit  multivariate-regression  model  is: 

PPR  =  10  £Eepk)0'728  (o;f°‘°39  Eq.  4.3 

F.  Shock  Tests 

1 .  Sand 

The  results  of  the  shock  tests  on  the  saturated  sand  are  given  in 
Chapter  III  in  Figures  3.4  and  in  Appendix  A.  The  best  fit  regression 
model  is  given  by  Equation  3.7  in  terms  of  peak  strain,  effective  stress 
and  relative  density.  The  best  fit  multivariate  regression  model  in  terms 
of  void  ratio  effective  stress  and  peak  strain  is: 

PPR  =  16  {Tepk)°'33  (o;f0-31  (Drf°-18  Eq.  4.4 

PPR  =  18  <Eepk)°'55  (c;f°-35  (Dr)'°-18  Eq.  4.4 

2.  Sand-Clayey  Silt  Mixture 

The  results  of  the  shock  tests  on  the  saturated  sand-clayey  silt 
mixture  are  given  in  Tables  4.7  and  Figures  4.10  and  4.11.  The  best  fit 
multivariate-regression  model  is  the  same  as  for  clayey  silt,  shown  in 
Equation  4.5.  These  two  soils  were  analyzed  together  since  there  were  no 
dramatic  differences  in  behavior. 
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3.  Clayey  Silt 

The  results  of  the  shock  tests  on  the  saturated  clayey  silt  at  a 
void  ratio  of  0.95  are  given  in  Table  4.7  and  Figures  4.10  and  4.11.  The 
best  fit  multivariate  regression  model  for  both  the  sand-clayey  silt 
mixture  and  the  clayey  silt  is: 

PPR  =  10.2  (£epk)0-728  Kf0*039  Eq.  4.3 


6.  Summary 

The  tests  results  indicate  that  the  residual  porewater  pressure 
increase  for  Monterey  No.  0/30  sand  is  approximately  the  same  for  both 
quasi-static  and  shock  loading.  The  residual  porewater  pressure  increase 
for  both  the  sand-silt  mixture  and  the  clayey  silt  is  greater  for  quasi¬ 
static  than  for  shock  loading. 
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Figure  4.2  Change  in  peak  porewater  pressure  as  a  function  of  the  change 
in  confining  pressure  in  terms  of  soil  type,  void  ratio,  and 
initial  effective  stress  for  all  static  tests. 


Table  4.2  List  of  Symbols  for  Figures  4.2  to  4.11 


Static  Tests 

Q  Monterey  No.  0/30  sand.  0R  =  8.61%.  effective  stress  =  172  kPa 
C3  Monterey  No.  0/30  sand,  0R  =  9.11%,  effective  stress  =  517  kPa 
A  Monterey  No.  0/20  sand.  DR  =  63.3?%.  effective  stress  = 


Monterey  No.  0/30  sand.  0R  -  63.80%,  effective  stress  5 
517  V.P& 

•*.  9Uty  sand,  e  -  0.669.  effective  stress  =  172  kPa 
+  Silty  sand,  e  =  0.669.  effective  stress  '  517  kPa 
®  Bonny  silt,  e  =  C.932,  effective  stress  =  172  kPa 
^  6o';ny  s»lt.  e  -  0-939.  effective  stress  =  517  kPa 


Oynam ic  Tests 

*  S;Uy  sand,  e  -  0.627,  effective  stress  =  172  kPa 
$  ijnry  silt,  e  =  0.935.  effective  stress  =  l?2  kPa 
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Figure  4.7 


Comparison  of  porewater  pressure  ratio  for  static  tests  as 
predicted  by  experimental  analysis  and  Veyera  (1935)  equation 
for  dynamic  tests  (Monterey  No.  0/30  sand)  utilizing  0D. 


SILT,  SAND  -  SILT  MIXTURE 


figure  4.10  Porewater  pressure  ratio  as  a  function  of  the  sum  of  the  peak 
sample  strains  in  terms  of  soil  type  and  void  ratio  for 
dynamic  tests. 
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Figure  4.11  Porewater  pressure  ratio  as  a  function  of  the  sum  of  the  peak 
sample  strains  in  terms  of  soil  type  and  void  ratio  for 
dynamic  tests. 
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Table  4.1.  Physical  Properties  of  Monterey  No.  0/30  Sand,  Bonny  Silt 
and  a  50-50  Mixture  of  These  Soils 


Monterey  No.  0/30  Bonny  Silt  50-50  Mixture 
Sand  Sand-Silt 

uses 


Classification 

SP 

ML 

SM 

Specific  Gravity 

2.65 

2.63 

2.64 

Particle  Size  Data 

°10 

0.29  Iran 

<  0.001  mm 

0.01  mm 

D30 

0.38  mm 

0.025  nro 

0.05  mm 

°50 

0.45  mm 

0.04  mm 

0.15  mm 

Cu 

1.65 

- 

34.1 

Cc 

1.00 

- 

0.55 

%  passing  #200  sieve 

0.06 % 

83.9% 

45% 

%  clay  size 

0% 

14 

6% 

%  sand 

99.94% 

16.1% 

5.5% 

%  silt 

0.05% 

69.9% 

39% 

LL 

0.06% 

25 

- 

PL 

- 

21 

- 

PI 

Relative  Density 

1700  kg/m3 

4 

- 

max  density 

— 

minimum  dry  density 

1470  kg/m3 

- 

- 

Proctor  Test 

Dry  Density 

“ 

1730  kg/m3 
@  14.8% 
water  content 

-wijw 
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Table  4.7  Dynamic  Test  Results  -  50%  Sand  -  50%  Silt  Mixture 
and  Bonny  Silt 


test  ID 

Effective 

Stress 

Void 

Ratio 

Impact 

upk 

(KPa) 

^pk 

(%) 

'  PPR 

50%  sand 

!72kPa 

.627 

1 

4298 

.08599 

.441 

50%  silt 

2 

4162 

.08328 

.666 

3 

23* 

.00046* 

.536 

4 

3438 

.06879 

.737 

5 

4342 

.08689 

.666 

6 

2986 

.05974 

.682 

Bonny 

i?2kPa 

.935 

1 

2126 

.05159 

.385 

Silt 

2 

1832 

.04445 

.441 

3 

1606 

.03896 

.544 

4 

3121 

.07573 

.666 

5 

2397 

.05816 

.678 

6 

2533 

.06146 

.668 

*  All  data  not  read  by  Biomation 
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V.  FIELD  EXPLOSIVE  TESTS  ON  PLACED  SAND 

This  chapter  summarizes  our  results  of  explosive  induced  planar 
and  spherical  stress  wave  tests  on  water  saturated  sand  conducted  by  Bretz 
(1988),  Hassen  (1988),  and  Schure  (1988). 

A.  Introduction 

We  conducted  a  series  of  field  explosive  tests  (21  spherical  and  6 
planar  detonations)  on  large  samples  of  placed  Poudre  Valley  sand.  The 
tests  were  designed  to  systematically  evaluate  porewater  pressure 
increases  and  liquefaction  of  water  saturated  cohesionless  soil  subjected 
to  explosive  induced  compressive  stress  waves.  The  soil  was  tested  at 
three  relative  densities  of  1,  50,  and  89  percent  for  the  spherical 
explosive  detonations  and  at  a  relative  density  of  89  percent  for  the 
planar  explosive  detonations.  The  transient  porewater  pressure,  total 
stress  and  acceleration  were  measured  with  transducers  and  recorded  with 
transient  data  recorders  during  the  passage  of  the  stress  wave  and  the 
residual  porewater  pressure  after  the  passage  of  the  stress  wave  was 
measured  with  a  piezometer.  Predictions  based  on  AFOSR  laboratory 
liquefaction  research  conducted  by  Veyera  (1985)  and  Hubert  (1986)  as 
described  in  Chapter  III,  as  well  as  the  empirical  models  given  in  Chapter 
II,  are  analyzed  and  compared  with  the  field  results. 

B.  Test  Site 

We  performed  the  field  explosive  tests  at  the  explosive  test  site 
located  at  the  Colorado  St?te  University  Engineering  Research  Center.  The 
facility,  developed  with  AFOSR  funding  has  a  State  of  Colorado  blasting 
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permit  and  allows  control  of  soil  conditions.  As  shown  in  Figure  5.1,  the 
large  soil  sample  is  located  below  the  regional  groundwater  table. 
Locating  the  sample  below  the  water  table  minimized  potential  stress  wave 
reflections  from  the  sample's  boundaries. 

C.  Instrumentation 

We  instrumented  the  sand  during  the  sand  placement  with  porewater 
pressure  transducers,  accelerometers-,  and  stress  gages  to  record  the 
transient  response  as  a  function  of  distance  and  charge  weight. 
Piezometers  were  utilized  to  obtain  the  late  time  porewater  pressure 
response.  Digital  transient  data  recorders  were  utilized  to  record  the 
response  of  the  embedded  instruments  and  high  speed  video  cameras  recorded 
the  tests  and  the  late  time  porewater  pressure  response.  Strain  gages 
were  used  to  determine  the  changes  in  soil  density.  The  exact  location  of 
the  instrumentation  is  given  by  8retz  (1988),  Hassen  (1988)  and  Schure 
(1988). 

D.  Soil  Properties 

The  tests  were  conducted  on  Poudre  Valley  sand,  which  is 
commercially  produced  by  crushing  gravel  obtained  from  the  Poudre  River. 
This  sand  has  an  angular  grain  shape,  a  Dgg  grain  size  of  0.52  mm,  a 

specific  gravity  of  2.68  and  is  classified  as  an  SP  in  the  USCS 
classification  system.  The  physical  and  index  properties  are  given  in 
Table  5.1  and  the  gradation  is  given  in  Figure  5.2.  Details  of  this 
material  are  given  in  Chapter  III,  Section  H.  Quasi-static  stress-strain 
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curves  for  the  sand  under  uniaxial  (kQ)  loading  conditions  are  shown  in 
Figures  5.3  and  5.4. 

E.  Soil  Placement  and  Saturation 

The  sample  size  was  4.4  m  in  diameter  by  1.7  m  high.  The  soil 

3 

volume  was  26  m  and  the  dry  mass  of  the  sample  was  approximately  45,000 
kg.  We  placed  the  sand  loose  in  0.3  meter  thick  lifts.  Each  lift  was 
compacted  (if  required)  with  a  vibratory  compactor  to  the  required  dry 
density.  After  the  final  layer  was  placed,  carbon  dioxide  (CO2)  was 

introduced  to  displace  the  air  in  the  soil  voids  and  to  aid  in  saturating 
the  soil.  After  the  air  was  displaced,  we  displaced  the  CO^  with  water  by 

upward  flushing  the  soil  with  warm  deaired  water.  Flushing  the  soil  was 
continued  for  about  two  weeks  until  seismic  velocity  indicated  full 
saturation  (Vc  s  1,500  m/sec). 

F.  Test  Procedure 

Once  the  sand  was  placed  and  saturated,  we  added  1.8  m  of  water 
above  the  top  of  the  sample.  Explosives  were  detonated  in  the  water  1.2  m 
above  the  sand  and  0.6  m  below  the  water  surface.  Seven  spherical  and  six 
planar  detonations  were  performed  on  the  saturated  sand  at  a  relative 
density  of  approximately  89  percent,  seven  spherical  detonations  at  a 
relative  density  of  approximately  1  percent  and  seven  spherical 
detonations  at  a  relative  density  of  approximately  50  percent. 
Piezoresistive  porewater  pressure  transducers  (ENDEVC0  Model  8511a-5kMl), 
resistive  total  stress  cells,  and  piezoresistive  accelerometers  were 


-84- 


util  i  zed  to  measure  the  soil's  porewater  pressure,  total  stress  and 
acceleration  during  the  passage  of  the  explosively  induced  stress  wave. 
Piezometers  were  used  to  measure  the  residual  porewater  pressure  response. 
After  each  explosive  detonation,  the  settlement  of  the  sample  was  measured 
and  checked  against  measurements  recorded  by  the  strain  gages.  After  the 
tests  were  completed,  the  water  was  pumped  out  and  the  final  location  of 
the  gages  and  other  objects  recorded. 

We  obtained  the  explosives  from  Buckley  Power  Company  of 

o 

Englewood,  Colorado.  Detonating  cord  (Primacord  manufactured  by  the 
Ensign-Bickford  Company)  was  used  for  the  planar  detonations.  A  water  gel 

D 

(Tovex  800  manufactured  by  Dupont),  was  used  for  the  smaller  spherical 

R 

detonations.  Tovex  800  has  an  energy  rating  of  894  calories  per  gram. 

G.  Planar  Stress  Wave  Test  Results 

Hassen  (1988)  conducted  explosive  field  tests  on  saturated  Poudre 
Valley  sand  which  had  been  placed  at  a  relative  density  of  90  percent.  To 
produce  a  planar  compressive  stress  wave,  a  7  m  diameter  grid  of 

D 

detonating  cord  (Primacord  )  was  placed  in  the  water  1.2  m  above  the 
sand's  surface.  Because  the  burn  speed  of  detonating  cord  is 
approximately  6,700  m/sec,  the  center  of  the  grid  was  lifted  to  ensure  a 
true  plane  stress  wave.  The  grid  was  center  detonated  with  an 
instantaneous  electric  blasting  cap.  The  detonating  cord  had  spacings  of 
0.3  m  or  0.6  m  depending  on  the  detonating  cord  explosive  rating  and 
explosive  density  required  to  produc  the  required  peak  stress. 


Total  explosive  mass  for  the  six  grids  detonated  ranged  from  0.531 

2 

to  2.55  kg,  giving  a  charge  density  ranging  from  0.15  to  0.73  gm/m  . 
Figure  5.5  shows  a  typical  porewater  pressure  time  history  as  the  stress 
wave  passes  through  the  sample  as  measured  by  the  pi ezoresi stive  porewater 
pressure  transducers.  Figure  5.6  shows  a  typical  porewater  pressure  ratio 
versus  time  as  measured  by  the  piezometers.  Table  5.2  and  Figures  5.7  and 
5.8  show  the  measured  peak  porewater  pressure,  upk,  and  the  peak  pore 

pressure  ratio,  PPR,  respectively,  as  a  function  of  scaled  distance.  The 
best  fit  equations  for  this  test  data  are: 

upk  =  30073  (-j^-)'0*342  Eq.  5.1 

PPR  =  1.11  (-4—  f0*0844  Eq.  5.2 

where  upk  is  the  peak  porewater  pressure  in  kPa,  R  is  distance  from  charge 

2 

in  meters,  and  wi  is  the  charge  intensity  in  gm/m  .  The  best  fit  equation 
for  the  pore  pressure  ratio  as  a  function  of  the  peak  compressive  strain 
is: 

PPR  =  1.48  Upk)°’29  Eq.  5.3 

where  epk  is  the  peak  strain  in  percent  calculated  from  Equations  2.7  and 

2.8. 


Spherical  Stress  Wave  Test  Results 


Bretz  (1988)  conducted  explosive  field  tests  on  saturated  Poudre 
Valley  sand  which  had  been  placed  at  a  relative  density  of  89  percent.  To 
produce  a  spherical  compressive  stress  wave,  point  charges  were  placed  in 


the  water  1.2  m  above  the  center  of  the  soil  sample  and  0.6  m  below  the 
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water  surface.  Six  charge  masses  ranging  from  0.03  kg  to  7.02  kg  were 
detonated.  Schure  (1088)  conducted  similar  tests  on  saturated  Poudre 
Valley  sand  on  loose  sand  (0R  -  1%)  and  medium-dense  sand  (DR  *  50%). 

Tables  5.3  and  5.4  and  Figures  5.9,  5.10  and  5.11  show  the  measured 

1/3 

PPR  as  a  function  of  scaled  distance  (R/W  )  and  peak  compressive  strain, 
epk*  ^e  kes*  fit  equation  for  the  peak  porewater  pressure  for  all 


spherical  tests  is: 


V  .  5°,°°o  (jjTTT)'1'5  »« 


Eq.  5.4 


The  best  fit  equations  obtained  by  Sretz  (1988)  for  the  soil  at  a  relative 
density  of  89  percent  are: 


PPR  =  1.15  (-JTJ-)"0'25 


Eq.  5.5 


PPR  =  1.23  (ep)c) 


0.167 


Eq.  5.6 

The  best  fit  equations  for  the  soil  at  a  relative  density  of  about  50 
percent  are: 

d  .i  oi 

Eq.  5.7 


PPR  =  57.2  (-173-f1’87 

w 


PPR  =  57.1  (epk) 


1.25 


Eq.  5.8 


The  best  fit  equations  obtained  by  Schure  (1988)  for  the  soil  at  a 
relative  density  of  about  1  percent  are: 


PPR  =  177  (— J73“>~2* 3 
W 


Eq.  5.9 


1.53 


PPR  =  174  (Cpk) 


Eq.  5.10 
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Figures  5.12,  5.13  and  5.14  present  the  field  data  for  spherical 
detonations  including  Equations  5.6,  5.8  and  5.10  for  the  pore  pressure 
ratio  as  a  function  of  peak  strain  at  relative  densities  of  89,  50  and  1 
percent,  respectively.  As  expected,  the  soils  at  higher  densities  require 
greater  peak  strains  to  liquefy. 


Sutwnan 


For  the  tests  at  a  relative  density  of  89  percent  conducted  by 


Hassen  (1988)  and  Bretz  (1988),  the  best  fit  equations  and  the  field  data 
for  the  residual  pore  pressure  increase  induced  by  explosive  induced 
planar  and  spherical  stress  waves  and  for  laboratory  shock  induced  planar 


stress  wave  are  shown  in  Figure  2.  Also  plotted  on  Figure  5.2  is  Equation 
3.7  presented  by  Veyera  {1985}  for  laboratory  plane  comoressive  shock 
loading  of  Monterey  No.  0/30  sand.  The  prediction,  based  on  the 
laboratory  derived  Equation  3.7,  over  predicts  the  pore  pressure  ratio  for 
both  tne  planar  (Eq.  5.3)  and  spherical  (Eq.  5.6)  explosive  detonations. 
The  spherical  explosive  detonations  result  in  a  higher  pore  pressure  ratio 
that  the  planar  detonations  producing  a  given  peak  compressive  strain. 
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Figure  5.1  Cross-section  of  explosive  test  pit  and  soil  sample. 


*4  tia  «n  HI  Kl  trim  «i»  U.s.  STD.  SIEVES 


Figure  5.2  Poudre  Valley  sand  gradation  curve. 
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Figure 


5.3  Stress-strain  curves  for  uniaxial  loading  of  Poudre  Valley 
sand  (arithmetic  scale). 

raUOK!  VALLEY  SAND 
1-0  LOAD— UNLOAD  CYCLE  ..  OR  —  BO  X 


5.4  Stress-strain  curves  for  uniaxial  loading  of  Poudre  Valley 
sand  (logarithmic  scale). 


Figure  5.10  Pore  pressure  ratio  for  spherical  shots  as  a  function  of 
scaled  distance  (0D  =  50%;  Schure,  1988). 
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Table  5.1  Indej  Yoperties  of  the  Poudre  Valley  Sand 


Unified  Soil  Classification  S. stem  Group  Symbol  SP 

(angular  grain  shape) 

Mean  Specific  Gravity,  Gs  2.68 

Particle  Size  Distribution  Data 
D1q,  mm  0.18 

O^q,  mm  0.25 

05g,  mm  0.52 

Cu,  (1)  3.85 

Cc,  (2)  0.95 


Relative  Density  Test  Data 

Dry  Unit  Weight,  kg/m3 
Maximum 
Minimum 

Void  Ratio 

Maximum  0.76 

Minimum  0.438 


Quartz  Mean  Specific  Gravity  2.65 


Quartz  Bulk  Modulus,  MPa 


1860 

1522 


cu  *  (D30)2/(D60  x  D10> 

<2>  Cc  =  W°10 


30680 
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Table  5.2.  Peak  Residual  Pore  Pressure  Increase  for  Poudre  Valley 
Sand  as  a  Function  of  Charge  Weight,  Distance  and  Soil  Density 
for  Planer  Tests  (Hassen,  1988} 


Test 

No. 

Total 

Charge 

Wt  (kg) 

Average  Charge 
Weight  per 

Area  (kq/m2) 

Peak  Rise  of 
Water  in 

Piezometer  (m) 

PPR 

Explosive 

Used 

Dense 

(relative 

density  approximately  89  percent) 

PI 

0.53 

0.15 

1.37 

0.86 

Primacord** 

3.8  g/m 

(61x61  cm  grid) 

P2 

0.74 

0.21 

1.55 

0.97 

D 

Primacord 

5.3  g/m 

(61x61  cm  grid) 

P3 

1.46 

0.42 

1.68 

1.05 

Primacord** 

10.7  g/m 
(61x61  cm  grid) 

P4 

1.46 

0.42 

1.39 

0.87 

Primacord** 

10.7  g/m 
(61x61  cm  grid) 

P5 

2.55 

0.73 

D 

Primacord 

10.7  g/m 
(30x30  era  grid) 

P6 

1.28 

0.26 

Primacord** 

5.3  g/m 

(30x30  cm  grid) 

Charge  Location  In  water  1.22  m  above  tank  and  0.61  m  below  water  surface. 


Piezometer  is  located  at  a  depth  of  1.67  m  below  soil  surface  and  at  a 
distance,  R,  from  the  charge  grid  of  2.89  m.  The  rise  of  water  required  for 
PPR  =  1  is  1.6  m. 


mb 
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Table  5.3.  Peak  residual  pore  pressure  increase  for  Poudre  Valley  sand  as  a 
("unction  of  charge  weight,  distance  and  soil  density  for  spherical  tests 

(Bretz,  1988). 


Test 

No. 


Charge  Peak  Rise  of  PPR 
Wt.  Water  in 

W  (kg)  Piezometer  (m) _ 


Scaled  Explosive 

Distance,  R/W1//3  Used 
(m/kg1/3) _ 


Dense  (Relative  Density  Approximately  89  Percent) 


SI 

0.03 

1.22 

0.65 

9.3 

Primacord 

S2 

0.08 

1.26 

0.69 

6.7 

Primacord 

S3 

0.25 

1.56 

0.86 

4.6 

Tovex  800 

S4 

0.76 

1.56 

0.84 

3.2 

Tovex  800 

S5 

2.25 

1.30 

0.70 

2.2 

Tovex  800 

S6 

0.76 

1.55 

0.83 

3.2 

Tovex  800 

S7 

7.02 

1.10 

0.65 

1.5 

Tovex  800 

Piezometer  is  located  at  a  depth  1.67  m  below  soil  surface  and  at  a 
distance,  R,  from  the  charge  of  2.89  m.  The  rise  of  water 
required  for  PPR  =  1  is  1.6  m. 


Charge  Location  in  water  1.22  m  above  tank  and  0.61  m  below  water  surface 
(S7,  1.22  m  above  tank  and  1  m  below  water) 
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Table  5.4  Peak  residual  pore  pressure  Increase  for  Poudre  Valley  sand  as  a 
function  of  charge  weight,  distance  and  soil  density  for  spherical  tests 
_ _ (Schure,  1988). 
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VI.  IN— SITU  EXPLOSIVE  TESTS  ON  ALLUVIAL  SAND 

This  chapter  summarizes  our  results  of  buried  point  source 
explosive  detonations  (spherical  stress  wave  tests)  on  an  alluvial  sand 
deposit  conducted  by  Jacobs  (1988). 

A.  Introduction 

We  conducted  a  series  of  in-situ  field  explosive  tests  on  an 
alluvial  sand  deposit  to  systematically  evaluate  explosive- induced 
porewater  pressure  increases  and  liquefaction  of  a  water  saturated  natural 
sand  deposit.  The  transient  and  long-term  porewater  pressure  response  was 
measured  and  compared  to  values  predicted  by  various  empirical  models 
developed  through  previous  laboratory  and  field  research  given  in  Chapter 
II. 

B.  Test  Site  and  Soil  Properties 

The  field  test  site  is  a  natural  braided  river  deposit  located  in 
the  South  Platte  River  approximately  4  km  north  of  Kersey,  Colorado.  We 
conducted  the  tests  on  a  sand  island  located  in  the  river  channel.  Two 
borings  and  two  cone  penetrometer  tests  indicated  that  the  first  3.5  m  are 
composed  of  a  dense  poorly  graded  sand  (SP)  having  an  in-situ  relative 
density  of  approximately  85  percent.  Below  3.5  m  is  silt  (ML)  with  a 
liquid  limit  of  40,  plastic  limit  of  25  and  a  natural  water  content  of  28 
percent.  Seismic  tests  indicate  that  the  soil  is  saturated  below  a  depth 
of  0.25  m  (wave  velocity  of  1,582  to  1,595  m/sec).  Shear  wave  velocities 
ranged  from  162  m/sec  at  a  depth  of  1  meter  to  277  m/sec  at  a  depth  of  3.5 
meters.  Shale  bedrock  is  known  to  be  located  at  a  depth  of  approximately 
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33  m  In  this  area.  Based  on  the  soil  samples  collected  at  0.61  m 
intervals  in  the  two  borings  on  the  sand  island,  the  soil  . igs  in  Figure 
6.3  were  developed  for  each  boring.  Laboratory  tests  were  performed  on 
the  recovered  samples.  Table  6.1  presents  the  results  for  the  upper  3.65 
m  sand  layer  and  Table  6.2  presents  the  results  for  the  underlying  silt 
layer.  Figure  6.1  shows  the  profile  at  the  island  based  on  the  two 
borings,  two  CPT  tests  and  seismic  data.  For  further  details  of  the  site, 
see  Rwebyogo  (1987)  and  Jacobs  (1988). 

C.  Instrumentation 

We  instrumented  the  site  at  the  locations  shown  in  Figures  6.2  and 
6.3  with  porewater  pressure  and  particle  velocity  transducers  to  record 
the  transient  response  as  a  function  of  distance  and  charge  weight. 
Piezometers  were  utilized  to  obtain  late  time  porewater  pressure  response. 
Digital  transient  data  recorders  were  utilized  to  record  the  transient 
response  and  high  speed  video  cameras  recorded  the  tests  and  the  excess 
late  time  residual  porewater  pressure  response.  Surveying  equipment  was 
used  to  determine  ground  surface  elevations. 

D.  Explosives  and  Test  Procedures 

The  explosives  we  utilized  were  obtained  from  Buckley  Power 

D 

Company  of  Englewood,  Colorado.  Detonating  cord  (Primacord  ,  manufactured 
by  the  Ensign-Sickford  Company)  was  used  for  the  smaller  detonations  and  a 

p 

water  gel  (Tovex  800  manufactured  by  Dupont)  was  used  for  the  larger 

D 

detonations.  The  energy  rating  of  Tovex  SCO  is  894  calories  per  gram. 


which  is  about  10  percent  less  than  TNT.  Instantaneous  electric  blasting 
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caps  were  used  to  initiate  all  detonations.  The  blasting  cap  utilized  had 
a  maximum  delay  of  2.0  micro  seconds  after  the  electric  energy  was  applied 
by  the  blasting  machine. 

Six  detonations  ranging  from  0.0045  to  9.06  kg  of  explosives  were 
detonated  at  a  depth  of  3.0  to  3.5  m.  Detonations  1  to  5  occurred  in 
borehole  number  one  and  the  final  detonations  occurred  in  borehole  number 
two. 

E.  Test  Results 

Table  6.3  lists  the  test  sequence,  explosive  type,  weight  and 
depth  of  explosive  detonated  during  the  test  program.  The  average 
compression  wave  velocity  from  the  explosive  to  the  water  pressure  and 
velocity  gages  was  1,680  m/sec. 

i.  Peak  Transient  Porewater  Pressure 

The  transient  porewater  pressure  was  measured  at  a  depth  of  3.5  m 
at  distances  of  6.1  and  12.2  m  for  detonations  1  to  5  and  at  15.5  and  19.2 
m  for  detonation  number  C.  Figure  6.4  shows  a  typical  porewater  pressure 
time  history  for  an  explosive  detonation  recorded.  Using  the  "cube  root" 
1/3 

scaling  law,  R/W  '  ,  the  peak  porewater  pressure  values  obtained  for  each 
detonation  at  each  porewater  pressure  transducer  are  plotted  on  Figure 
6.5.  The  equation  representing  the  statistically  best  fit  line  using  a 
linear  least  squares  analysis  for  the  alluvial  site  studied  is 

_  r„  „„„  ,  R  ,-2.38 

V  *  50-093  y? r> 


Eq.  6.1 
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where  u  ,  is  the  peak  porewater  pressure  In  kPa,  R  is  the  distance  from 
pk 

the  detonation  in  m,  and  W  is  the  charge  mass  in  kg.  Figure  6.6  compares 
the  recorded  data  and  Equation  6.1  to  other  equations  given  in  the 
literature  given  in  Chapter  II. 

1i.  Peak  Particle  Velocity 

The  particle  velocity,  monitored  at  the  groundwater  table,  was 
measured  in  the  transverse,  longitudinal  and  vertical  directions.  The 
maximum  vector  sum  of  the  three  directions  is  presented  as  the  peak 
particle  velocity.  The  data  is  plotted  in  Figure  6.7.  The  statistically 
best  fit  line  for  the  alluvial  site  is 


Vak  ■  8'745  <^7T>'2-06 


Eq.  6.2 


where  Vpeajc  is  the  peak  particle  velocity  in  m/sec.  Figure  6.8  compares 

the  recorded  data  and  Equation  6.2  to  other  equations  in  the  literature 
given  in  Chapter  II. 


iii.  Settlement  Survey 

The  survey  results  revealed  that  there  was  no  significant 
settlement  in  the  area  surrounding  the  explosive  detonations  for 
detonations  1  to  5.  Detonation  number  6  (9.1  kg)  produced  a  ground  heave 
of  5  cm.  These  observations  are  consistant  for  a  dense  sand  deposit. 


iv.  Residual  Porewater  Pressure  and  Porewater  Pressure  Ratio 
Residual  porewater  pressure  was  measured  by  five  piezometers.  The 
elevation  rise  in  each  standpipe  was  recorded  by  using  video  recorders. 
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Tile  results  are  presented  as  a  dimensionless  number  known  as  the  porewater 
pressure  ratio,  PPR,  defined  in  Chapter  2.  For  the  depth  of  the 
piezometer  opening  (3.5  m),  the  effective  stress  for  the  groundwater  table 
elevation  at  the  time  of  the  tests  was  approximately  38  kPa.  This  means 
that  a  rise  of  387  cm  in  any  of  the  piezometer  stand  pipes  would  yield  a 
pore  pressure  ratio  of  1.0,  defining  a  condition  of  liquefaction.  The  PPR 
time  histories  for  each  piezometer  for  the  sixth  detonation  (9.10  kg)  is 
given  in  Figure  6.9.  The  time  histories  for  the  other  detonations  are 
given  in  Jacobs  (1988).  The  piezometers  have  a  response  time  (natural 
frequency)  of  approximately  one  second. 

The  measured  peak  residual  porewater  pressure  ratio,  PPR,  for  all 
detonations  is  plotted  in  Figure  6.10.  The  statistically  best  fit 
equation  for  the  measured  data  is 


PPR  =  2.60  (-TTrf1-366 


Eq.  6.3 


Because  of  the  high  permeability  of  the  alluvial  sand  at  the  site,  the 
residual  porewater  pressure  dissipated  very  fast  (less  than  five  seconds). 
The  piezometer  system,  with  its  response  time  of  approximately  one  second, 
was  unable  to  measure  the  peak  residual  porewater  pressure  immediately 
following  the  passage  of  the  stress  wave.  Utilizing  Terzaghi's 
consolidation  theory,  the  actual  peak  residual  porewater  pressure  is 

?PIW)  =  4-8  f^TT)"1'48  E"'  6  4 


where  PPR(pea^)  is  the  residual  porewater  pressure  ratio  just  after  the 

passage  of  the  stress  wave  and  before  consolidation  occurs.  Equation  6.4 
is  plotted  in  Figure  6.11.  This  equation  predicts  liquefaction  (PPR  =  1) 

at  a  scaled  distance,  of  2.9  m/kg*^.  Lyakhov  (1961)  equation  in 
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Chapter  II  predicted  liquefaction  at  a  scaled  distance  between  2  and  8, 
depending  on  depth  of  the  explosive  soil  type  and  soil  density.  Ivanov 
(1967)  equation  predicts  liquefaction  at  a  scaled  distance  of  less  than  6 
1/3 

m/kg  .  Studer  and  Kok's  (1980)  predicts  liquefaction  at  a  scaled 

distance  of  2.8  m/kg*^.  Equation  6.4  fits  Studer  and  Kok's  (1980) 
equation  the  best. 

1/3 

Utilizing  a  scaled  distance  of  2.9  m/kg  '  in  Equation  6.4,  the 
peak  particle  velocity  required  to  cause  liquefaction  at  the  South  Platte 
River  Site  is  0.98  m/sec.  Dividing  the  peak  particle  velocity  by  ere 
measured  compressive  wave  velocity  of  1,680  m/sec  leads  to  a  peak 
compressive  strain  of  0.058%  required  to  cause  liquefaction  at  the  site. 
Solving  Veyera's  (1985)  equation  given  in  Chapter  II  at  a  relative  density 
of  89%  and  an  effective  stress  of  38  kPa  leads  to  the  following 

PPR  =  2.38  (epk)°'331  Eq.  6.5 

where  epk  is  the  peak  compressive  strain  in  percent.  As  shown  in  Figure 

6.12,  at  a  strain  of  0.06,  a  PPR  of  0.93  is  predicted  by  Veyera's  equation 
which  is  very  close  to  the  peak  strain  causing  liquefaction  at  the  site. 

F.  Summary 

The  measured  peak  transient  porewater  pressure  induced  by  an 
explosive  detonation  in  the  South  Platte  alluvial  soil  deposit  was  found 
to  be  similar  magnitude  to  the  predicted  by  Drake  and  Little  (1983)  and 
Drake  and  Ingram  (1981)  empirical  equations.  The  empirical  equations 
developed  by  Ivanov  (1957),  Lyakhov  (1961)  and  Cole  (1948)  over-predict 
the  peak  transient  porewater  pressure  in  the  South  Platte  River  deposit. 
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The  measured  peak  particle  velocities  Induced  by  explosive- 
detonations  in  the  South  Platte  River  deposit  was  found  to  be  slightly 
greater  in  magnitude  than  predicted  by  empirical  equations  presented  by 
Drake  and  Little  (1983)  and  Drake  and  Ingram  (1981). 

The  calculated  PPR,  utilizing  the  permeability  and  drainage 
conditions  at  the  South  Platte  River  deposit  compared  well  to  the 
empirical  equations  by  Studer  and  Kok  (1980)  based  on  scaled  distance  and 
to  Veyera's  (1985)  equation  based  on  explosive  induced  strain. 


SOUTH  PLftTTE  RIVER  SOIL  PROFILE  8-9 
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INSTfcUriEKTftTION  LAYOUT  FOR  DETONATIONS  i-S 
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SP1  -  SP5  =  settlement  plates  one  through  five 
PI  -  P5  =  piezometers  one  through  five 
PV1  -  PV2  «  particle  velocity  gauges  one  and  two 
PPT1  -  PPT2  =  porewater  pressure  transducers  one  and  two 


Figure  6.2  Instrumentation  layout,  plan  view,  for  detonations  one  through 
five  which  utilizes  borehole  one  as  the  explosive  detonation 
point. 
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INSTRUtlENTaTION  LAYOUT  FOR  DETONATION  6 
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SP1  -  SP5  =  settlement  plates  one  through  five 
PI  -  P5  =  piezometers  one  through  five 
PV1  -  PV2  =  particle  velocity  gauges  one  and  two 
PPT1  -  PPT2  =  porewater  pressure  transducers  one  and  two 


Figure  6.3  Instrumentation  layout,  plan  view,  for  detonations  six  which 
utilizes  borehole  two  as  the  explosive  detonation  point. 
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Figure  6.7  Peak  particle  velocity  data  and  the  empirical  equation 
representing  the  line  of  best  fit  through  the  field  data. 
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Figure  6.8  Peak  particle  velocity  data  for  the  South  Platte  River  filed 
site,  the  best  statistically  best  fit  line  through  the  data 
and  other  equations  given  in  the  literature. 
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Figure  6.9  Residual  porewater  pressure  expressed  as  the  dimensionless 
porewater  pressure  ratio,  (PPR  =  u/a'},  vs.  time  after 
detonation  number  six,  {9.1030  kg). 
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Figure  6. 10  Peak  residual  porewater  pressure  ratio 'for  all  detonations  and 
the  empirical  equation  representing  the  line  of  best  fit 
through  the  data. 
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Figure  6.11  Calculated  peak  residual  porewater  pressure  ratio  at  t  =  o+  for 
all  detonations  and  Equation  6.4. 
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figure  6.12  Statistically  best  fit  equations  for  peak  residual  porewater 
pressure  ratio  predictions  based  on  peak  strain. 


Table  6.1  Summarized  laboratory  test  results  for  the  upper 

3.65  m  sand  layer 


Soil  Classification  (USCS,  ASTM  02487  and  0442) 

SP  -  rounded  to  sub-rounded  coarse  gravely  poorly  graded  sand. 
%  passing  #200  <  1% 

Bgo  "  2  nw 


Minimum  and  Maximum  Oensity  (ASTM  04254  and  04253) 


Ymax 

e  • 
rmn 


=  77  KN/m3 
=  0.373 


Ymin 

emax 


a  16.48  KN/m3 
=  0.564 


Tahir  6.2  Summarized  laboratory  test  results  for  the  silt  layer, 
intercepted  at  a  depth  of  3.65  m 

Soil  Classification  (USCS,  ASTM  D2487  and  D422) 

ML  -  low  plasticity  silt 
Atterberg  Limits  (ASTM  04318) 

Liquid  Limit  -  40 
Plastic  Limit  -  25 
Water  content  -  28.5% 

Plasticity  Index  -  15 
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Table  6.3  Test  sequence,  charge  type,  weight  and  depth  of  explosive  for 
the  six  detonations  fired  during  the  test  program. 


Test 

number 

Date 

July.  1987 

Depth  of  charge 
placement 
(meters) 

Type  of 
Explosive 

Charge 

mass 

(kilograms) 

1 

14th 

3.5 

blasting  cap 

0.0045 

2 

15th 

3.5 

Primacord 

0.0453 

3 

15th 

3.0 

Tovex 

0.091 

4 

15th 

3.0 

Tovex 

0.453 

5 

16th 

2.8 

Tovex 

2.264 

6 

16th 

3.0-3. 5 

Tovex 

9.103 
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VII.  ANALYSIS  OF  SHOCK  AND  EXPLOSIVE  INDUCED  LIQUEFACTION 


This  chapter  presents  our  empirical,  analytical  and  theoretical 
models  to  predict  liquefaction. 


a.  Empirical 

As  given  in  Chapters  II  and  III,  several  empirical  equations  have 
been  proposed  to  predict  explosive  induced  porewater  pressure  increases 
and  liquefaction.  The  major  ones  are  presented  below: 

PPR  «  1.65  +  0.64  In  (W1/3/R)  Eq.  2.1 

PPR  =  16.3  (2epk}°'331  (o^f0,308  {DR)~‘179  Eq.  2.9 

R  *  2.5  W1/3  for  PPR  =  1  Eq.  2.2 

max 

R  =  fc,  W1/3  for  PPR  *  1  Eq.  2.3 

ITiaX  X 

where  PPR  is  the  pore  pressure  ratio  defined  in  Equation  3.4,  W  is  the 

charge  mass  in  kg,  R  is  the  radial  distance  from  the  charge  in  m,  Ep^  is 

the  peak  compressive  strain  in  percent,  is  the  effective  stress  in  kPa 

and  0R  is  relative  density  in  percent.  Our  field  explosive  tests  on 


placed  sand  (Chapter  V)  and  our  field  explosive  tests  on  in-situ  sand 
(Chapter  VI)  result  in  the  following  empirical  equations: 


PPR  =  30073  (4-) 
PPR  -  1.15  (^3) 

PPR  =  57.2  (-*7,) 


-0.342 


-0.25 


-1.87 


planar  for  Dr=89% 
spherical  for  Dt>=89% 

spherical  for  Dr=50% 


Eq.  5.1 
Eq.  5.5 

Eq.  5.7 
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PPR  =  117  (-f73)“2‘3  spherical  for  Dr=l%  Eq.  5.9 

PPR  =  4.8  {R/W173)*1-48  spherical  for  Dr=90%  Eq.  6.4 

PPR  =  2.38  Upk)°’331  spherical  for  0^=90%  Eq.  6.5 
8.  Analytical 

Analytical  methods  to  model  explosive  induced  porewater  pressure 
increases  and  liquefaction  in  the  two-phase  water  saturated  sands  have 
been  developed  by  Veyera  (1985),  Charlie  et  al .  (1987),  and  Awad  (1988). 
These  models  utilize  Biot's  (1956,  1962)  theory  solved  numerically  by 

either  the  finite  element  method  or  the  finite  difference  method.  Because 
of  the  steep  wave  front  and  the  high  velocity  of  stress  wave  propagation 
in  water  saturated  soil,  very  small  time  steps  and  small  element  size  are 
required  for  numerical  stability,  convergence  and  accuracy  (Zienkiewicz 
and  Bettess,  1982).  For  one-dimensional  modeling,  this  does  not  present 
major  computational  problems  but  two-  and  three-dimensional  modeling 
requires  very  large  computer  memory  and  considerable  computer  time.  For 
example,  Awad's  (1988)  axi-symmetric  finite  element  model  takes  over  60 

D 

minutes  to  run  on  a  C0C/2Q5  super  computer. 

Although  these  analytical  models  have  had  only  limited 
verification,  these  analytical  models  have  predicted  explosive  induced 
porewater  pressure  increases  and  liquefaction  of  water  saturated 
cohesionless  soils. 

C.  Theoretical 

Bretz  (1988)  developed  a  closed  form  solution  to  theoretically 
predict  radial  and  tangential  strain  and  stress  due  to  the  detonation  of 
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spherieal  charges  in  water  saturated  soil.  The  method  models  water 
saturated  soil  as  an  elastic  water-solid  mixture  having  zero  shear 
strength  and  assuming  that  the  soil  particle  contact  can  be  ignored. 
Assuming  that  the  soil  skeleton  follows  the  same  strain  path  as  the 
mixture,  linear  and  non-linear  stress-strain  characteristics  of  the  soil 
skeleton  can  be  followed  both  during  the  loading  and  unloading  phase.  The 
method  closely  matches  the  results  of  field  explosive  tests  on  saturated 
sands  reported  by  Bretz  (1988).  The  basic  equations  utilized  by  Bretz 
(1988)  are: 

(1  +  e)  8  B 

Bm  =  — g — B -  (Wood  Equation)  Eq.  8.1 

°S  °W 

where  is  the  bulk  modulus  of  the  mixture;  8s  and  8w  are  bulk  moduli  of 
the  granules  and  water,  respectively;  and  e  is  the  void  ratio. 

°r  *  (1  VM)(1  --2M)  CU  '  v)  §7  +  2u  7]  Eq.  8.2 

°r  *  (1  +  lv  f+  rJ  8'3 

where  or  and  are  radial  and  tangential  stresses,  u  is  Poisson's  ratio; 

u  is  radial  displacement;  r  is  radius  at  which  u  occurs;  and  E  is  modulus 
of  elasticity. 


c 


v 


Eq.  8.4 


where  cv  is  the  coefficient  of  consolidation;  ug  is  excess  porewater 

pressure;  z  is  a  coordinate  in  the  verticle  direction,  and  t  is  time. 
This  equation  is  the  consolidation  equation,  and  it  was  used  with  a  finite 
difference  solution  to  calculate  excess  porewater  pressure  after  passage 
of  the  dynamic  stress  wave,  i.e.,  in  the  t  +  1  regime.  Since  soils  are 
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multiphase  materials,  the  soil  skeleton  is  inelastic  and  excess  porewater 
pressure  dissipation  occurs  with  time,  a  true  closed  form  theoretical 
model  for  multiphase  soil  would  be  very  complex  and  may  be  impossible  to 
formulate.  Theoretical  modeling  of  the  soil  particle  contacts  to  better 
understand  the  nonlinearity  of  granular  material  is  needed. 
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VIII.  SUMMARY 

Our  field,  laboratory  and  theoretical  research  indicate  that  the 
destruction  potential  of  an  explosion  may  be  greatly  magnified  if 
detonated  in  water  saturated  granular  soils.  While  blast-induced 
liquefaction  may  not  necessarily  damage  a  facility  structurally,  it  may 
render  it  unusable.  Empirical  models  are  given  that  can  be  used  to 
estimate  liquefaction  potential  as  a  function  of  density,  effective  stress 
and  applied  compressive  strain.  One  of  the  models  uses  an  empirical 
scaling  law  for  explosive  loadings  to  predict  the  extent  of  porewater 
pressure  increases  in  the  field  from  buried,  contained  charges  in 
saturated  soils.  A  numerical  analysis  that  considers  the  saturated  soil 
as  a  two-phase  medium  is  presented.  The  analysis  accounts  for  the 
nonlinear,  inelastic  behavior  of  the  soil  skeleton  and  has  shown  that 
liquefaction  is  dependent  upon  the  constrained  modulus  of  the  soil 
skeleton.  Results  agree  with  the  experimental  observations  of  peak  and 
long-term  porewater  pressure  responses. 

The  results  of  our  study  indicate  the  following. 

1.  Liquefaction  can  be  induced  by  single  and  multiple  blasts. 

2.  Liquefaction  can  be  induced  at  distances  much  greater  than 
those  associated  with  structural  damage. 

3.  Long  term  increases  in  residual  porewater  pressures  can  be 
induced  by  compressive  shock  wave  loadings  when  the  peak 
particle  velocity  exceeds  0.075  m/s. 

4.  Liquefaction  can  be  induced  in  loose  saturated  sand  by  a 
single  compressive  shock,  wave  when  the  peak  particle  velocity 
exceeds  0.75  m/s. 
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5.  Soils  at  higher  initial  effective  stress  and  higher  initial 
relative  density  require  more  energy  to  produce  liquefaction. 

6.  Destruction  potential  of  an  explosive  charge  may  be  greatly 
mags, 'fled  if  detonated  in  water-saturated  soils. 

7.  Liquefaction  occurs  because  of  compressive  strain  induced  by 
the  compression  stress  wave,  but  liquefaction  occurs  after  the 
stress  wave  passes. 

8.  Liquefaction  occurs  because  loading-unloading  of  the  porewater 
is  elastic  but  the  soil  skeleton  is  not. 

An  explosive  detonated  in  a  soil  having  a  high  liquefaction 
potential  could  cause  damage  disproportionate  to  the  energy  released. 
Documented  occurrence  of  blast-induced  liquefactior  is  available  in  the 
open  literature.  Although  considerable  work  remains  to  be  done  in 
projecting  this  information  into  a  comprehensive  method  of  predicting 
liquefaction  for  actual  or  hypothetical  blasts,  the  data  indicate  that 
residual  porewater  pressure  increases  should  not  occur  in  soils  subject  to 
strains  of  less  than  O.OC5  percent.  Transient  and  quasi -static  tests 
indicate  that  residual  porewater  pressure  increases  and  liquefaction  are 
not  strain-rate  sensitive  for  sand  but  are  strain-rate  sensitive  for  silt. 
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APPENDIX  A.l.  MONTEREY  NO.  0/30  QUARTZ  BEACH  SAND 

1.  Skeleton  Stress-Strain  Curves 

Static  one-dimensional,  confined,  compression  tests  (Hendron, 
1963;  Whitman  et  al.,  1964)  were  performed  on  air  dry  samples  of  Monterey 
No.  0/30  sand.  The  results  of  tests  were  used  to  obtain  stress-strain 
information  for  the  soil  skeleton  and  to  determine  the  constrained  modulus 
of  the  skeleton  to  be  used  in  calculating  the  theoretical  porewater 
pressure  response  (C-parameter)  and  for  modeling  of  the  residual  porewater 
pressure.  To  simulate  the  initial  stress  conditions  in  the  experimental 
investigation,  a  corresponding  initial  effective  stress  was  applied  to 
each  sample.  The  compressive  strain  values  were  referenced  to  the  initial 

applied  stress.  Each  sample  was  loaded  and  unloaded  in  increments  two 

times  to  develop  the  stress-strain  relationship  for  the  soil  skeleton. 
From  these  results,  a  constrained  modulus  for  loading  and  unloading  was 
determined. 

Samples  were  tested  at  relative  densities  of  40  and  80  percent. 
Each  sample  was  tested  at  initial  effective  stresses  of  86  kPa  and  690 

kPa.  The  skeleton  stress-strain  curves  for  Monterey  No.  0/30  sand  are 

shown  in  Figures  A. 2  through  A. 5.  The  results  show  that  the  soil  skeleton 
stiffness  increases  with  increasing  initial  effective  stress  and  density. 
Hysteresis  between  the  loading  and  unloading  curves  decreases  with 
increasing  initial  effective  stress  and  density. 

2.  Static  C-Parameter  Response 

8efore  loading  of  each  sample,  the  porewater  pressure  response  was 
checked  to  determine  the  degree  of  saturation.  This  was  done  by 
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increasing  the  confining  pressure  on  the  sample  and  monitoring  the 
sample's  porewater  pressure  response  without  drainage.  The  ratio  of  the 
sample  porewater  pressure  response  to  the  increase  in  confining  pressure 
is  termed  the  “C-parameter1*  (Lamfae  and  Whitman,  1969)  for  a  one¬ 
dimensional  confined,  compressive  loading  of  a  saturated  soil  with 
undrained  conditions.  A  ratio  of  one  indicates  a  saturated  sample  and 
values  less  than  one  indicates  that  the  sample  is  not  saturated  or  has  a 
stiff  soil  skeleton.  For  an  initial  effective  stress  of  86  kPa,  a  ratio 
of  one  was  consistently  obtained.  However,  a  C-parameter  of  less  than  one 
was  obtained  for  higher  initial  effective  stresses.  Since  the  preparation 
and  saturation  process  was  identical  for  each  sample,  it  was  assumed  that 
the  porewater  pressure  ratios  obtained  indicated  of  a  saturated  condition. 
An  examination  of  the  compressive  stress  wave  propagation  velocity  through 
samples  verified  this  assumption.  The  measured  compressive  stress  wave 
velocities  were  close  to  1,500  meters  per  second  in  all  samples 
investigated.  This  is  the  value  that  would  be  expected  for  saturated 
conditions. 

Throughout  the  experimental  investigation  it  was  noted  that  the 
porewater  pressure  ratio  varied  in  a  predictable  manner  with  variations  in 
effective  stress  and  relative  density  (Figure  A. 8).  In  considering  this 
observation  and  those  previously  discussed,  it  is  believed  that  the 
porewater  pressure  ratio  response  noted  can  be  attributed  to  changes  in 
the  soil  skeleton  stiffness  which  increases  with  increasing  initial 
relative  density  and  effective  stress.  Accordingly,  all  samples  were 
considered  to  be  saturated. 
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3.  Pressure-Time  Histories 

The  pressure-time  histories  represent  the  porewater  pressure 
transducer  responses  to  applied  shock  loadings  as  function  of  time.  They 
include  both  the  peak  and  long-term  response  for  the  confining  pressure 
and  the  sample  porewater  pressure.  Selected  pressure-time  histories, 
representative  of  the  behavior  observed  in  this  experimental 
investigation,  are  shown  in  Figures  A. 7  and  A. 8.  A  summary  of  numerical 
results  from  the  pressure-time  histories  of  all  samples  investigated  is 
given  in  Tables  A. 2  through  A. 12. 

Figure  A. 7  shows  the  experimental  result  for  the  relative 

density  series  at  an  effective  stress  of  172  kPa.  Figure  A. 8  shows  the 
experimental  result  for  the  "80%"  relative  density  series  at  an  effective 
stress  of  172  kPa.  The  confining  pressure  and  sample  porewater  pressure 
responses  have  been  plotted  together  on  each  figure.  The  "series" 
designation  for  relative  density  has  been  used  to  group  together  data 
having  approximately  the  same  r.ative  density.  The  designations  include 
data  that  is  within  10  percent  greater  than  the  series  number  (including 
the  series  number).  For  example,  a  "40%"  series  designation  would  include 
all  data  for  a  relative  density  from  40  percent  to  49  percent. 

The  pressure-time  histories  at  i  indicative  of  the  system  response  ! 

during  and  after  loading.  On  each  ft.  re,  the  traces  of  the  confining 
pressure  and  sample  porewater  pressure  follow  each  other  closely  in  their 
response  trends.  The  two  curves  are  slightly  offset  from  one  another  in  ! 

the  time  domain  due  to  the  relative  locations  of  each  transducer.  The 
sample  peak  porewater  pressure  is  greater  than  the  applied  stress  peak 
values  for  each  impact.  ! 


In  all  cases,  the  confining  pressure  transducer  response  returned 
to  its  original  baseline  value  once  the  compressive  stress  wave  energy  had 
dissipated.  The  confining  pressure  should  return  to  its  original  value  if 
the  system  is  not  allowed  to  drain.  The  residual  excess  porewater 
pressure  indicated  by  the  sample  transducer,  was  above  its  original 
baseline  value  after  each  loading  and  continued  to  increase  with  each 
successive  impact.  The  response  of  the  sample  porewater  transducer  was  as 
expected  since  an  increase  in  the  residual  excess  porewater  pressure 
should  be  maintained  for  undrained  conditions.  Liquefaction  occurs  when 
the  residual  excess  porewater  equals  the  effective  stress  which  is  also 
when  the  back  pressure  plus  the  residual  excess  porewater  pressure  equals 
the  confining  pressure. 

The  information  from  the  pressure-time  history  records  was  used 
for  the  analysis  of  data  prssented  in  Tables  A. 2  through  A. 12.  The 
porewater  pressure  response  was  evaluated  as  a  function  of  initial 
effective  stress,  initial  sample  density,  peak  compressive  strain. 
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APPENDIX  A. 2.  TYNDALL  QUARTZ  BEACH  SAND 

1.  Static  Compression  Tests  on  the  Soil  Skeleton 
Static  compression  tests  were  not  run  on  these  samples. 

2.  Static  C-Parameter  Response 

The  saturation  process  was  easily  accomplished  and  the  C-parameter 
was  .92,  .88  and  .77  for  relative  densities  of  55,  63  and  73  percent. 
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APPENDIX  A. 3.  POUDRE  VALLEY  GRANITIC  SAND 

1.  Skeleton  Stress-Strain  Curves 

Samples  of  air  dry  soil  were  tested  at  initial  relative  densities 
of  20,  40,  60  and  80  percent.  The  skeleton  stress-strain  curves  are  shown 
in  Chapter  V. 

2.  Static  C-Parameter  Response 

The  C  parameter  varied  between  0.95  for  loose  samples  under  a  low 
effective  stress  to  0.5  for  dense  samples  under  a  high  effective  stress. 
Compression  wave  velocity  through  the  sample  of  approximately  1,500  m/sec 
confirmed  that  the  sample  was  saturated. 

3.  Shock  Induced  Settlement 

In  Figure  A. 11,  the  settlement  ratio  is  plotted  versus  the  initial 


relative  density. 
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APPENDIX  A 4 .  POUORE  VALLEY  GRANITIC  FINE  SAND  AND  GRAVEL 


1.  Skeleton  Stress-Strain  Curves 

Static  compression  tests  were  not  run  on  these  samples. 

2.  Static  C-Parameter  Response 

The  static  C-parameter  ranged  from  0.85  to  0.81  for  the  sand  and 
from  0.98  to  0.89  for  the  gravel. 


APPENDIX  A. 5.  ENIWETOK  CORAL  BEACH  SAND 


One-dimensional  static  compression  tests  were  performed  on  dry 
samples  of  Eniwetok  coral  sand  to  evaluate  the  stress-strain  behavior  of 
the  soil  skeleton  and  to  estimate  its  crushing  potential.  The  samples,  at 
40  and  80  percent  relative  density  with  an  initial  low  and  high  confining 
pressure  (150  and  517  kPa)  were  loaded  and  unloaded  incrementally  twice  to 
generate  the  stress-strain  relationship  of  the  soil  skeleton.  Figures 
A. 15  to  A. 20  show  the  variation  in  stress-strain  for  each  test.  The  low 
density,  low  effective  stress  test  is  more  compressible  than  the  high 
density,  high  effective  stress  test. 

In  Figures  A. 15  and  A. 16,  for  low  effective  stresses,  crushing  of 
the  particles  can  be  recognized  by  the  reverse  curvature  of  the  stress- 
strain  curve  (Lambe  and  Whitman,  1969).  The  result  of  these  coral  sand 
compression  tests  show  a  definite  concave  trend  toward  the  strain  axis  on 
the  loading  path.  Popping  sounds  could  be  heard  during  loading  indicating 
also  a  moderate  amount  of  particle  fracturing. 

2.  Static  C-Parameter  Response 

For  all  tests,  the  C-parameter  value  was  less  than  one  (Figure 
A. 21).  For  low  effective  stress  and  low  relative  density  tests,  the 
porewater  pressure  response  was  usually  higher  than  for  the  tests  with 
denser  sands  and  higher  effective  stresses.  This  difference  may  be 
attributed  to  the  changes  in  soil  skeleton  stiffness  since  the 
compressibility  of  a  soil  skeleton  will  decrease  with  increases  in 
effective  stress  and  in  relative  density  (Lee  et  al . ,  1969).  Veyera 
(1985)  reported  similar  porewater  pressure  response  was  recorded  for 
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saturated  dense  sands.  Accordingly,  most  tests  were  considered  saturated. 
Secondary  porosity  of  the  individual  sand  grain,  where  the  air  may  not 
have  been  totally  removed  and  system  compliance  may  also  explain  the  lower 
C-parameter.  Measured  compression  wave  velocities  of  1500  meters  per 
second  or  greater  through  the  sample  indicate  full  saturation. 

3.  Shock  Induced  Settlement 

Shock  induced  settlement  as  a  function  of  relative  density  is 
given  in  Figure  A. 17 


Figure  A.l  Grain  size  distribution  for  Monterey  No.  0/30  sand  (Muzzy, 1983; 
Chari ie  et  al . ,1984) . 
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Figure  A. 2  Skeleton  stress-strain  curve  for  Monterey  No.  0/30  sand  at  Dr 
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Figure  A. 3  Skeleton  stress-strain  curve  for  Monterey  No.  0/30  sand  at  D 


40%  and  o'  =  690  kPa. 
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Figure  A. 10  Grain  size  distribution  of  Poudre  Valley  sand. 
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Figure  A. II  Settlement  ratio  for  Poudre  Valley  Sand. 
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Figure  A. 17  Skeleton  stress-strain  curve  for  Eniwetok  coral  sand  at  D; 


Figure  A. 18  Skeleton  stress-strai n  curve  for  Eniwetok  coral  sand  at  D 
80%  and  o'  =  517  kPa. 
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Figure  A. 21  Static  C-Parameter  response  as  a  function  of  the  effective 
stress  for  Eniwetok  coral  sand. 
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Figure  A. 22  Settlement  ratio  for  each  test  versus  initial  relative 
density. 


Table  A.l  Physical  Properties  of  Monterey  Mo.  0/30  Sand 
(Muzzy,  1983;  Charlie  et  al.,  1984) 
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Table  A. 2  Peak  Porewater  Pressures  and  Peak  Compressive  Strains  for 
Monterey  No.  0/30  Sand  for  Low  Impact  Stress  Loading  and  Dr  =  “0%"  Series 
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Table  A. 3  Peak.  Porewater  Pressures  and  Peak  Compressive  Strains  for 
Monterey  No.  0/30  Sand  for  Low  Impact  Stress  Loading  and  Df  =  "20%“  Series 
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Table  A. 4  Peak  Porewater  Pressures  and  Peak  Compressive  Strains  for 
Monterey  No.  0/30  Sand  for  Low  Impact  Stress  Loading  and  Dr  =  “40%  “  Series 
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Table  A. 5  Peak  Porewater  Pressures  and  Peak  Compressive  Strains  for 
Monterey  No.  0/30  Sand  for  Low  Impact  Stress  Loading  and  =  “60%u  Series 

(First  Set) 
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Table  A. 6  Peak  Porewater  Pressures  and  Peak  Compressive  Strains  for 
Monterey  No.  0/30  Sand  for  Low  Impact  Stress  Loading  and  Dr  =  "60%“  Series 

(Second  Set) 
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Table  A. 7  Peak  Porewater  Pressures  and  Peak  Compressive  Strains  for 
Monterey  No.  0/30  Sand  for  Low  Impact  Stress  Loading  and  Dr  -  "80%  "  Series 
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Table  A. 8  Peak  Porewater  Pressures  and  Peak  Compressive  Strains  for 
Monterey  No.  0/30  Sand  for  High  Impact  Stress  Loading  and  =  “0%“  Series 
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Table  A. 9  Peak  Porewater  Pressures  and  Peak  Compressive  Strains  for 

Monterey  No.  0/30  Sand  for  High  Impact  Stress  Loading  and 

D„  =  “20%H  Series 
r 


Tott  I.P. 

(CP*) 

0 

t 

(*> 

Impott 

UP*) 

<*) 

fP*<5 

*28X1 

172 

21.* 

1 

3*17 

.01332 

1.015 

2 

4510 

.0*995 

.994 

3 

5033 

.109*4 

.994 

4 

47J0 

.10349 

.9*4 

£20X2 

S4S 

27.5 

1 

5400 

.1171* 

.994 

2 

5453 

.12271 

.994 

3 

574* 

.12517 

.124 

4 

3421 

.07425 

.437 

St  0X4 

(M 

22.1 

1 

57<t 

.12*24 

1.035 

2 

<274 

.11733 

1.013 

3 

<247 

.13*73 

.994 

4 

500 

.1113* 

.9*4 

Mot*:  (1)  ImiikI 

<2)  C»le*l»t*d 

-  Urn  lo*4i»j»  *«r<  4qm  tot  • '  *  14 


Table  A. 10  Peak  Porewater  Pressures  and  Peak  Compressive  Strains  for 

Monterey  No.  0/30  Sand  for  High  Impact  Stress  Loading  and 

0  =  “40%"  Series 
r 


T«st  1.0. 

UP*) 

D 

r 

<« 

I*p**t 

.  <« 
*4 

UP*) 

P* 

(*) 

PPX<1> 

140X1 

172 

44.2 

1 

4524 

.09550 

.504 

2 

<447 

.14032 

.»* 

3 

49C2 

.14571 

.*29 

4 

4419 

.13974 

.171 

340X2 

545 

44.7 

1 

44*3 

.0*444 

.134 

2 

3474 

.07727 

•  919 

3 

4015 

.01440 

.919 

4 

4249 

.0*974 

.921 

*40X4 

<90 

44.7 

1 

3*45 

.09013 

.717 

2 

5570 

.11701 

.911 

3 

497* 

.10440 

.932 

4 

7211 

.15140 

1.074 

Not*:  (1)  ■**»«*< 

-  No  1 0*41*14  ttn  4*a*  tot  *  14  IP* 
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Table  A. 11  Peak  Porewater  Pressures  and  Peak  Compressive  Strains  for 

Monterey  No.  0/30  Sand  for  High  Impact  Stress  Loading  and 

0  =  "60%“  Series 
r 


Tait  1.8. 


•  'a  B(  iMftat 

CEP»>  ft) 


pk 

UP*) 


m 


u> 


scon 

1T2 

<1.1 

1 

4156 

.01504 

.834 

2 

5211 

.11684 

1.166 

3 

3422 

.07000 

1.166 

4 

4326 

.08351 

1.143 

860X2 

245 

44.1 

1 

<022 

.12265 

.662 

2 

1117 

.10422 

.745 

3 

5594 

.11401 

.909 

4 

«1* 

.12259 

.994 

*60X4 

(to 

<2.5 

1 

53*4 

.12245 

.624 

2 

56$1 

.11607 

.704 

3 

5771 

.11790 

.787 

4 

5599 

•11438 

.787 

Hat*:  <1>  lunni 

(1)  C»U*Ut»4 


-  Ha  mi  tom  tot  •  H  Di 


Table  A. 12  Peak  Porewater  Pressures  and  Peak  Compressive  Strains  for 

Monterey  No.  0/30  Sand  for  High  Impact  Stress  Loading  and 

D  =  “8Q%"  Series 
r 


T*»t  X.D.  o’e  9^  Xoptct 

UP«>  ft) 


UP.) 


(2) 


ft) 


m. 


u> 


560X1 

172 

81.5 

1 

4867 

.13525 

1.154 

2 

41 IJ 

.08241 

1.204 

3 

5824 

.11470 

1.196 

4 

(923 

.13(47 

1.230 

54012 

345 

83.3 

i 

3191 

.16076 

.332 

2 

6913 

.13(37 

.497 

3 

7040 

.13805 

.364 

4 

8413 

.1(376 

.581 

wax  4 

<90 

83.8 

2 

5966 

.11(88 

.427 

3 

57  J7 

.11240 

.745 

3 

5717 

.11240 

.794 

4 

8371 

.16400 

.911 

Not*:  (1)  M*ti*r*4 

(21  C*)«*Ut«4 


N*  lot* i*i ■  **r*  4***  ter  *'  «  W  Oi 
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Table  A. 13  Experimental  Porewater  Pressure  Ratio  Predictor  Models 


fed*  i 


U)  PPK 


<»>  PP*  *  «< 


f(E  ...  (VW)-) 
<U.?«)  | - 1 

r  r  .  .  }•»’ 
*,J>  U'>  (Be)| 


.i7» 


(*>  m 


<9.411  CC  «  .)'3iS  (^l-238 


si. a 


61.X 


17. « 


.204 


.113 


•  207 


(d> 

m  - 

<12.71) 

“V314 

Dtr-M2 

43.* 

.1*1 

(«) 

PFX  • 

(11.33) 

<«’)“*305  (VMT14* 

9 

39.9 

.203 

II) 

PMt  - 

(14.30) 

«vm 

C»')"*308  (Dr,--»» 

63.* 

.1*7 

Note:  -  and  Dr  are  in  percent 

-  a1  is  in  kPa 

o 

-  D  =  relative  density 

-  VDP  =  Volume  Oecrease  Potential  =  e  -  e  . 

nn  h 

-  e  =  void  ratio  =  volume  of  voids/volume  of  solids 
2 

R  =  the  coefficient  of  determination 
S  =  the  standard  error  of  estimate 
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Table  A. 14  Tyndall  Florida  Beach  Sand  Characteristics  and  Tests  Performed 


Florida  Beach  Sand  Characteristics 


G  *  2.65  (quartz) 

B  {modulus)  *  30680  HPa  {quartz) 

W  1290  k3/n’ 


'sax 


(method  by  Kolbuszewski ,  1348) 
=  1724  kg/m 5 

{method  by  Kolbuszewski,  1948) 


Shock  Test  Parameters 

C-para«eter 
Cannon  pressure 


Tests  Performed 


Relative  Density 
Effective  stress 


s  .77  to  .32 
*  45.5  to  50.0  KPa 


*  50%.  60%.  70% 

=  345  KPa  {50  psi) 


USCS  Classification  SP  Particle  Size 

0l8  *  0.17  arc 
Dj„  1  0.21  m 
050  =  0.22  wra 
0,,,  =  0.31  ora 
Cu  =  1.47 
Cc  =  1.04 


Table  A. 15  Results  of  Tests  Performed  on  Tyndall  Beach  Sand 


* 


Test  ID 

e«‘ 

(kPa) 

Dr 

(%5 

Impact 

V 

(kPa) 

Epk 

{%) 

PPR 

50/50 

345 

54.7 

1 

995 

.0214 

.038 

f 

2 

973 

.0210 

.392 

3 

3438 

.0743 

.899 

4 

1945 

.0420 

.965 

5 

2375 

.0512 

.929 

6 

4478 

.0968 

.965 

60/50 

345 

62.6 

1 

1040 

.0232 

.526 

c 

2 

3822 

.0852 

.993 

3 

135 

.0030 

.953 

4 

4139 

.0923 

.992 

5 

4S8 

.0111 

1.022 

6 

4207 

.0938 

1.016 

70/50 

345 

72.8  . 

1 

68 

.0014 

.016 

2 

1674 

.0349 

.203 

• 

3 

498 

.0103 

.251 

4 

4455 

.0922 

.325 

5 

588 

>0122 

.732 

6 

1809 

,0377 

.834 

•I 


i 

j 

* 


Table  A. 16  Poudre  Valley  Sand  Characteristics  and  Tests  Performed 


Poudre  Valiev  Sand  Characteristics 


6s  *  2.68 

''min  *  •i49P  (93.0  f5cfi 

v  *  I860  kg/m1  (116.3  pcf) 


USCS  Classification  SP  Particle  Size 


Shock  Test  Parameters 

C-parameter  =  .5  to  .9 

Cannon  pressure  =  45.5  to  49.6  KPa  • 

Tests  Performed  (P  tests) 

Relative  density  (Dr)  =  20%,  40%,  60%,  80% 

Effective  stress  («(,*)  *  86  KPa  and  172  KPa 

(12.5  and  25  Psi) 


Diei  «  0.18  mm 

DI0  =  0.25  mm  One  test  at  80%  Df  and  517  KPa  (75Psi)  5 

O50  *  No  test  was  Perforw««1  for  *  relative  density  of  100% 

Cu°  -  3.85 

Cc  =  0.95 


Table  A. 17  Results  of  Tests  Performed  on  Saturated  Poudre  Valley  Sands 


Test  ID 

On* 

(KPa) 

°r 

<*) 

Impact 

Upk 

(KPa) 

tPk 

(V 

PPR 

1 

1193 

.025 

.798 

20/12.5 

85 

26.6 

2 

1220 

.025 

1.061 

3 

1827 

.033 

1.060 

4 

3523 

.075 

1.298 

1 

2699 

.055 

.670 

20/25 

172 

27. 

2 

3523 

.075 

.673 

3 

3730 

.079 

.691 

4 

1537 

.070 

.931 

1 

1241 

.025 

.788 

40/12.5 

86 

46. B 

2 

2261 

.045 

.924 

3 

1103 

.022 

1.002 

4  - 

2868 

.058 

1.062 

1 

136 

,002 

.108 

40/25 

172 

42.9  2 

1103 

.022 

.801 

3 

1220 

.025 

1-061 

4 

1289 

.326 

1.126 
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Table  A. 17  Continued 


Test  ID  o9*  Dr 

(KPa)  (%) 


60/12.5  86  66.5 


60/25  172  65.5 


80/12.5  86  81.1 


80/25  172  81.1 


80/75  517  80.9 


Impact 

Up* 

(KPa) 

Cpk 

(14) 

PPR 

1 

179 

.003 

.080 

2 

226 

.004 

.272 

3 

1310 

.024 

.878 

4 

1333 

.025 

1.042 

1 

2640 

.050 

.331 

2 

1310 

.024 

.353 

3 

1379 

.026 

.380 

4 

723 

.018 

.529 

1 

2420 

.043 

.789 

2 

2054 

.036 

.796 

3 

1537 

.027 

.796 

4 

1537 

.027 

.827 

1 

2889 

.051 

.536 

2 

2389 

.051 

.674 

3 

68 

.001 

.680 

4 

2641 

.047 

.801 

1 

2620 

.046 

.663 

2 

2551 

.045 

.766 

3 

2730 

.048 

.855 

4 

2620 

.046 

.915 

5 

2303 

.041 

.953 

6 

2847 

.050 

.955 

Table  A. 18  Physical  Properties  of  the  Poudre  Valley  Fine  Sand 

and  Gravel  Material 


Fine  usd 

Gravel 

uses  ication 

SP 

GP 

Spuciific  Gravity 

2. 68 

2.68 

Particle  Sire  Data: 

»10 

0*  10  ££2 

5.00  as 

d30 

0.16  sa 

6. 80  ms 

°50 

0.20  KB 

8. 90  ms 

B60 

0.23  KB 

12.00  &si 

1.43 

1.50 

ced) 

1.11 

1.07 

Relative  Denaity  Test  Data 

Drv  Unit  Weight: 

Maxiaua 

1655  kg/*3 

1676  Kg/si3 

Hinlaun 

14SE  kg/»3 

1527  Kg/*3 

Void  Ratio: 

Maxi run 

0.838 

0.757 

Miniaum 

0.618 

0.599 

e  A. 19  Peak  Porewater 

Pressures  and 

Peak  Compressive 

Strains 

for  the  Poudre  Valley  Fine  Sand, 

°r  = 

"30%"  Series 

Test  I.D  o0'  Dr 

Impact  Up* i 1 1 

fpk<2)  PPR(3) 

kPa  % 

kPa 

% 

C~para»eter  ■  0.85 

S133  207  37.4 

1 

68 

.00159 

.ill 

2 

203 

.00479 

.331 

3 

792 

.01864 

.781 

4 

588 

.01385 

.897 

5 

57 

.00054 

.914 

6 

1040 

.02450 

1.105 

C-paraceter  *  0.81 

S135  345  35.8 

1 

543 

.01282 

.398 

2 

316 

.00748 

.596 

181 

.00427 

.663 

Notes  (1)  Measured 

4 

159 

.00375 

.730 

(2}  Calculated 

5 

203 

.00481 

.796 

(3)  Calculated 

6 

272 

.00642 

.928 

167 


Table  A. 20  Peak  Porewater  Pressures  and  Peak  Compressive  Strains 
for  the  Poudre  Valley  Fine  Sand,  Dr  =  "70%"  Series 


fe»t  i. a  «„•  Dj 

UP«>  <»} 

In pact 

«pk 4 1  * 
(kf**} 

U) 

PPRO) 

C-par4»«C€r  -  0.82 

SI?  3  20?  77.3 

X 

452 

.00576 

•  333 

2 

475 

.03025 

» 539 

3 

475 

.01025 

.664 

4 

203 

.00439 

.672 

5 

1357 

.02929 

.995 

C-^nnw  -  0.  81 

6 

113 

.00244 

.  995 

sm  3*5  •n.j 

1 

701 

.01507 

.467 

2 

543 

.01167 

r  530 

2 

ioie 

.02168 

.  666 

4 

1016 

.02168 

.795 

5 

90S 

.02139 

.660 

6 

1176 

.02529 

.951 

Table  A. 21  Peak  Porewater  Pressures  and  Peak  Compressive  Strains 
for  the  Poudre  Valley  Gravel  Material,  =  "30%''  Series 


Test  I.D  oQ’  Dr 
(kPa)  (%) 


Impact 


upktn  4pltt2)  PFR<3) 
(kpa)  (%) 


C -parameter  «  . 98 
S233  207  39.7 


C-parameter  «  .94 
S235  345  38.4 


l 

757 

.01700 

.556 

2 

181 

.00411 

.652 

3 

339 

.00772 

.667 

4 

521 

.01185 

.819 

5 

973 

,  .00411 

1.043 

6 

950 

.02215 

1.108 

1 

1086 

.02481 

.662 

2 

724 

.01654 

.728 

3 

1131 

.02583 

.861 

4 

678 

.01550 

.927 

5 

1154 

.02635 

.993 

6 

701 

.01602 

.993 

I 


41 
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Table  A. 22  Peak  Porewater  Pressures  and  Peak  Compressive  Strains 
for  the  Poudre  Valley  Gravel  Material,  Dr  =  "70% "  Series 

Test  I.D  o0'  nr  I aspect  up)t(1*  -ipX{2)  PPR<3* 


(XFa)  (!) 

OtPa) 

{%) 

C-parameter  *  .93 

S273  207  78.6 

1 

498 

.01054 

.374 

2 

152 

.00309 

.452 

3 

639 

.01366 

.635 

4 

769 

.01628 

.721 

5 

1131 

.02394 

.992 

6 

860 

.01820 

1.028 

C-parameter  •  .89 

S275  345  78.4 

1 

1063 

.02237 

.513 

2 

113 

.00238 

.596 

3 

1040 

.02189 

.729 

4 

882 

.01855 

.788 

S 

1131 

.02379 

.862 

6 

1402 

.02950 

.994 

Table  A. 23  Physical  Properties  of  the  Eniwetok  Coral  Sand 


USCS  Classification 

SP 

Specific  Gravity 

2. SC 

Particle  Size 

Di. 

.32  mm 

DJ0 

.43  mm 

D,o 

.48  mm 

o„ 

.65  nsn 

Coefficient  of  Uniformity  Cy 

1.66 

Coefficient  of  Curvature  C£ 

1.09 

Relative  Density  (dry  unit  weight) 

Maximum 

1705  kg/m 

Minimum 

1509  kg/m 

Void  Ratio 

Minimus* 

.60S 

Maximum 

.316 

Bulk  Modulus  of  CaCOj  (Calcite)  * 

67500  MPa 

Bulk  Modulus  of  CaCOj  (Aragonite)** 

70000  MPa 

Note:  *  from  Goodman  (1976) 
**  from  Clark  (1966) 
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Table  A. 24  Peak  Porewater  Pressure  and  Peak  Compressive  Strain  for 
Eniwetok  Coral  Sand,  *  "0%"  Series 


Test  # 

°e‘ 

(KPa) 

D_  Impact 

r 

w 

v1 

(KPa) 

£p*! 

<%) 

PPft1 

S-00/25 

172 

6.1 

l 

22.75 

.0004 

.030 

2 

701.22 

.0182 

1.100 

3 

339.23 

.0077 

1.088 

4 

1017.70 

.0231 

1.070 

5 

565.39 

.0128 

1.063 

6 

1583.09 

.0358 

1.190 

S-00/50 

345 

10.0 

1 

0.00 

.0010 

.000 

2 

180.65 

.0040 

.203 

3 

113.08 

.0025 

.260 

. 

'  4 

339.23 

.0076 

.453 

5 

90.32 

.0020 

.457 

6 

791.54 

10179 

.676 

S-00/75 

517 

2.2 

1 

3008.29 

.0684 

1.015 

2 

2284.31 

.0519 

1.104 

3 

3460.60 

.0788 

1.104 

4 

3008.29 

.0685 

1.104 

5 

3777.08 

.0860 

1.104 

6 

2985.53 

.0680 

1.104 

Peak 

Porewater 

Pressure 

and 

Peak  Compressive 

Strain  for 

Eniwetok  Coral  Sand,  Dr  =  "2Q%"  Series 


Test  # 

(KPa) 

°r 

<%) 

Impact 

vl 

(KPa) 

Sk1 

(*) 

PPR1 

5-20/25 

172 

23.5 

1 

0.50 

.0001 

.001 

2 

1900.26 

.0422 

.977 

3 

1854,75 

.0411 

.947 

4 

1130.78 

.0251 

.954 

5 

1877.51 

.0416 

.923 

6 

950.13 

.0211 

.929 

S-2C/50 

345 

27.2 

1 

2329.82 

.0514 

.986 

2 

1379.69 

.0304 

.970 

3 

2329.82 

.0514 

.934 

4 

2239.50 

.0495 

.928 

„ 

•  5 

3483.35 

.0769 

.935 

6 

3SS3.42 

.0673 

.952 

5-20/75 

517 

26.2 

1 

2374.64 

.0524 

.851 

2 

1741.68 

.0385 

1.0CS 

3 

1176.29 

.0260 

1.101 

4 

1442.90 

.0539 

1.148 

5 

1153.53 

.0254 

1.144 

6 

1425.20 

.0315 

1.136 
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Table  A. 26  Peak  Porewater  Pressure  and  Peak  Compressive  Strain  for 

Eniwetok  Coral  Sand,  D  =  “40%“  Series 
_ _ _ _ ’  r  _ _ _ _ 


Test  1 

«»* 

Dr 

Impact 

V 

*pk! 

PPR’ 

(KPa) 

<%) 

{KPa) 

(%) 

S-40/25 

172 

47.8 

1 

542.64 

.0117 

.442 

2 

1312.12 

.0281 

.662 

3 

610.90 

.0132 

.794 

4 

814.30 

.0174 

.795 

5 

1221.10 

.0262 

.933 

6 

814.30 

.0174 

1.032 

S-40/50 

345 

39.7 

1 

1289.36 

.0280 

.596 

2 

2578.73 

.0559 

.927 

3 

3415.09 

.0740 

1.059 

4 

837.05 

.0181 

1.059 

5 

3799.83 

.0824 

1.062 

6 

950.13 

'.0206 

1.059 

$-40/75 

517 

43.8 

1 

316.48 

.0074 

.180 

2 

2691.80 

.0580 

.748 

3 

1017.70 

.0220 

.795 

4 

2736.62 

.0590 

.883 

5 

3211.69 

.0693 

.971 

6 

2352.57 

.0507 

.971 

Peak 

Porewater 

Pressure  and 

Peak  Compressive 

Strain  for 

Eniwetok  Coral  Sand 

•  Dr  = 

"60%”  Series 

Test  # 

®o' 

Dr 

Impact 

“dR1 

EpkJ 

PPR3 

(KPa) 

(*) 

(KPa) 

(%> 

S-60/25 

172 

59.7 

1 

2261.56 

.0476 

.929 

2 

1741.68 

.0367 

.929 

3 

994.95 

.0209 

.927 

4 

791.55 

.0167 

.920 

5 

2917.96 

.0614 

.929 

6 

2374.64 

.0499 

.928 

$-60/50 

345 

65.0 

1 

3528.86 

.0733 

.201 

2 

3302.01 

.0690 

.411 

3 

3935.67 

.0823 

.464 

* 

4 

4184.57 

.0875 

.502 

c 

4048.74 

.0846 

.465 

6 

3234.44 

.0676 

.464 

S-60/75 

517 

63.1 

1 

3437.85 

.0704 

.785 

2 

2872.46 

.0603 

.873 

3 

3257.20 

.0682 

.917 

4 

3550.92 

.0744 

.969 

5 

4184.57 

.0877 

1.002 

6 

2601.48 

.0545 

1.047 
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Table  A. 28  Peak  Porewater  Pressure  and  Peak  Compressive  Strain  for 


Eniwetok  Coral  Sand,  D„  =  "SOX'1  Series 

r 


Test  # 

o#‘ 

(KPa) 

°r 

{%) 

Impact 

u  J 
oli 

(KPa) 

epk* 

(%) 

PPRJ 

5*80/25 

172 

73.1 

1 

45.07 

.0010 

.259 

2 

1221.10 

.0252 

.511 

3 

2080.91 

.0429 

.673 

4 

3095.54 

.0640 

.882 

5 

7.96 

.0004 

.920 

6 

2239.50 

.0462 

.927 

S-80/50 

345 

75.0 

1 

3279.95 

.0675 

.463 

Z 

4546.56 

•093S 

.582 

■  '3 

90.24 

.0018 

.592 

4 

3890.15 

.0800 

.556 

5 

4478.30 

'  .0921 

.530 

6 

4365.22 

.0898 

.529 

S-80/75 

517 

78.3 

1 

3257.20 

.0666 

.839 

2 

1922.33 

.0394 

.923 

3 

3008.29 

.0615 

.971 

4 

3302.01 

.0676 

.971 

5 

2261.56 

.0463 

1.021 

6 

2307.07 

.0473 

1.044 

Table  A.  29  Peak  Porewater  Pressure  and  Peak  Compressive  Strain  for 


Eniwetok  Coral 

Sand, 

100%"  Series 

Test  *  o,1 

Dr 

Impact 

Upk‘ 

V* 

PPR3 

CKPa) 

(%) 

( KPa ) 

W 

S-100/25  172 

105.7 

1 

1990.59 

.0389 

.706 

2 

2261.56 

.0442 

1.266 

3 

2714.56 

.0530 

1.037 

4 

2397.39 

.0468 

.811 

5 

2442.90 

.0767 

.811 

6 

2374.64 

.0463 

.802 

$-100/50  345 

103.5 

1 

67.57 

.0013 

.067 

2 

859.80 

.6168 

.395 

•  3' 

4433.48 

.0869 

.664 

4 

1605.84 

.0315 

.724 

5 

1130.78 

.0221 

.730 

6 

2261.56 

.0441 

.775 

S— 1 00/7  5  517 

103.5 

1 

1130.78 

.0221 

-593 

2 

1877.51 

.0367 

.576 

3 

2050.91 

.0408 

.664 

4 

1425.20 

.0279 

.629 

5 

1153.53 

.0226 

.663 

6 

3279.95 

.0643 

.685 

Note:  1  Measured 


3  Calculated 


1  Calculated 
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B.l  RESEARCH  OBJECTIVES 

The  primary  objective  of  our  research  was  to  systematically 
evaluate  the  behavior  of  saturated  granular  soils  subjected  to  shock  and 
explosive  loadings.  Secondary  objectives  included  developing 
experimental,  empirical,  analytical  and  theoretical  methods  to  better 
understand  and  evaluate  blast  induced  liquefaction. 

B.2  SIGNIFICANT  ACCOMPLISHMENTS 

Our  development  of  laboratory  shock  and  field  explosive  facilities 
are  major  accomplishments.  These  facilities  have  the  capability  of  being 
upgraded  for  higher  stress  levels  and  for  conducting  different  types  of 
transient  tests  on  saturated  and  unsaturated  soils.  The  current  testing 
of  level  soil  deposits  can  be  expanded  to  sloping  ground,  soil -structure 
interaction,  and  pile  foundation  response.  The  experimental  testing  has 
allowed  us  to  develop  empirical  models  and  to  start  to  develop  and 
evaluate  theoretical  and  analytical  models.  Our  research  demonstrates 
that  the  destruction  potential  of  an  explosion  may  be  greatly  magnified  if 
detonated  in  water  saturated  granular  soils.  While  blast-induced 
liquefaction  may  not  necessarily  damage  a  facility  structurally,  it  may 
render  it  unusable.  Blast-induced  liquefaction  can  cause  late  time 
decreases  in  the  soil's  shear  strength  that  produces  damage 
disproportionate  to  the  amount  of  explosive  used  and  ground  motions 
inconsistent  with  previous  experience. 

The  results  of  our  study  indicate  the  following. 

1.  Liquefaction  can  be  induced  by  single  and  multiple  explosive 
induced  compressive  wave  loadings. 
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2.  liquefaction  can  be  induced  at  distances  from  explosions  much 
greater  than  those  associated  with  structural  damage. 

3.  Fairly  long  term  increases  in  residual  porewater  pressures  can 
be  induced  by  compressive  shock  wave  loadings  when  the  peak 
particle  velocity  exceeds  0.075  m/s  or  the  peak  strain  exceeds 
0.005  percent. 

4.  Liquefaction  can  be  induced  in  loose  saturated  sands  by  a 
single  compressive  shock  wave  when  the  peak  particle  velocity 
exceeds  0.75  m/s  or  the  peak  strain  exceeds  0.05  percent. 

5.  Soils  at  higher  initial  effective  stress  and  higher  initial 
relative  density  require  more  energy  to  produce  liquefaction. 

6.  Liquefaction  occurs  as  a  result  of  compressive  strain  induced 
by  the  compression  stress  wave,  but  liquefaction  occurs  during 
unloading  and  after  the  stress  wave  passes. 

7.  Liquefaction  occurs  because  loading-unloading  of  the  porewater 
is  elastic  and  reversible,  but  loading-unloading  of  the  soil 
skeleton  is  not  an  elastic-reversible  process. 
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8.7  CONSULTATIONS  TO  GOVERNMENT  AGENCIES 


ONA  (blast-induced  liquefaction) 

AFuil  (blast- Induced  liquefaction) 

AFESC  (blast- induced  liquefaction) 

GNR  (blast-induced  liquefaction) 

WES  (blasting  near  dams) 

Bureau  of  Reclamation  (blasting  near  dams) 
Bureau  of  Mines  (blasting  near  dams) 


S.8  NEW  DISCOVERIES  AND  INVENTIONS 


We  have  observed  and  recorded  blast-induced  liquefaction  in 
laboratory  shock  tests  and  field  explosive  tests.  Empirical,  analytical 
and  theoretical  methods  to  predict  and  better  understand  shock  induced 
porewater  pressure  have  been  developed  and  evaluated.  We  have  discovered 
that  peak  strain  controls  development  of  residual  porewater  pressure  and 
that  liquefaction  occurs  upon  unloading.  Loading  rate  is  not  important  in 
generation  of  residual  porewater  pressure  in  sands  but  is  important  for 
silty  soils. 

We  have  applied  for  a  patent  for  a  piezovane  developed  to  evaluate 
the  potential  for  blast-induced  liquefaction  leading  to  large  soil  strains 
(RTC  Disclosure  No.  044-0144-88) . 


